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Chapter 1 

Introduction 



The discovery of the first charmonium state [A"*" 74a I A"*" 74b] in November 1974 led to a revo- 



lution in particle physics, especially in hadron spectroscopy. It promoted the quark model and 
non-Abelian gauge field theories over competing scenarios such as the 'democracy of hadrons' 
|DRGG75] ■ leading to the prevailing picture of particle physics. The following years of research 
resulted in the subsequent observation of other charmonium resonances and the discovery of the 
first bottomonium states in 1977 [H"'"77 . Since that time a lot of theoretical effort has been di- 



rected toward understanding hadron spectra and properties within quark models. Quarkonium 
systems, i.e. bound states of a heavy quark and antiquark, have played a particularly important 
role in helping to establish the gauge field theory of quarks and gluons and their interactions, 
Quantum Chromodynamics (QCD). 

QCD features two remarkable properties. First, asymptotic freedom implies that at very high 
energies and momenta, quarks and gluons interact only weakly and act as quasifree particles 
|GW731 IPol73j . Second, confinement presumably results from the fact that at low energies the 
force between quarks increases with their distance, so that quarks are always tied into hadrons 
and cannot be removed individually. 

Confinement, a consequence of the non-Abelian structure of the theory, makes it hard to 
calculate quantities for bound states within QCD as one cannot apply perturbative QCD. By 
analogy with positronium, and given the large masses of the charm and bottom quarks, non- 
relativistic phenomenological potential models have been applied as tools for quarkonium spec- 
troscopy. To accommodate the properties of QCD these models, e.g. lEGK+781 IRic791 IBTSlj . 
are based on a short range part motivated by perturbative QCD and a phenomenological long 
range part accounting for confinement. 

After more than a decade of quarkonium physics, the discovery of new resonances above and 
below the DD-threshold, often referred to as the Charmonium Renaissance, has revived the 
interest in this field [B+04[ [Swa06] . 

Prom the first potential models for quarkonium in the late 70's to nowadays there have been 
many new developments in the field of quarkonium spectroscopy. Among others rigorous QCD 
evaluations |TY94l[TY95] . lattice simulations [O+021 IdP+Ol |JOO+06[ [KK06j or Effective Field 



Theory (EFT) calculations [ BBP9n[ IBPSV991 IBPSVOnl IBPSVOH fBSynil IBPSVOH] are applied 
to gain more fundamental insight. Nevertheless phenomenological potential models still play an 
important role in heavy meson spectroscopy. Unfortunately some of the more refined models of 
the mid 80's to nowadays tend to use potentials with a high degree of uncontrolled phenomeno- 
logical input IGRR821 085l IGSR891 IZVOR951 IEJ^'G03l IRROT] . 

In this context the task of this diploma thesis is to build up successively a potential model, 
which yields quantitative quarkonium spectra by using a minimum of parameters and minimal 
phenomenological input. Consequently the study concentrates on input with fundamental ori- 
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Chapter 1. Introduction 



gin, either already used, checking whether essential for the model or not, or so far unaccounted 
for. A second task of this work is to investigate the relation between radius and mass of the 
g^-states. 

This thesis is organized as follows: Chapter [2] gives an introduction to basic concepts of quarko- 
nium potential models. 

In Chapters [3] & H] two versions of the basic Coulomb-plus-Linear model, the first one in- 
cluding retardation corrections and the second one without them, are discussed. In addition the 
relation between mass and radius of qq-states within the model is investigated in Chapter [3l 

Chapter [5] introduces the concept of an induced t-channel interaction and discusses a model 
including this interaction. 

A model based on a perturbative QCD potential including processes up to order 0{a'^), and 
a linear confining potential is investigated in Chapter [6l 

In Chapter [7] the model arising from a combination of the induced interaction and a higher 
order perturbative QCD potential is analyzed and discussed in detail. Further the radius-mass 
relationship is investigated for this model. 

The thesis ends with a short summary in Chapter [5J In this chapter the results of the 
previous chapters are briefly summarized and a conclusion and outlook is given. 

The models discussed in this work are denoted with acronyms. The acronym reflects the principle 
structure of the model: 

• [CpL-B] 

Coulomb-plus-Linear model with perturbative QCD part due to Breit interaction; 

• [CpL-3] 

Coulomb-plus-Linear model with perturbative QCD part due to one-gluon exchange and 
momentum dependent terms according to option 3 in Appendix [D| 

• [CpLpL3] 

Coulomb-plus-Linear model with induced interaction, perturbative QCD part due to one- 
gluon exchange and momentum dependent terms according to option 3 in Appendix [Dj 

• [CpL*-2] 

Coulomb-plus-Linear model with perturbative QCD part to order 0{a'^) and momentum 
dependent terms according to option 2 in Appendix \D\ 

• [CpL*pI-2] 

Coulomb-plus-Linear model with induced interaction, perturbative QCD part to order 
O(a^) and momentum dependent terms according to option 2 in Appendix [Pl 

As common in this field of physics all equations are presented in rationalized units with 

h = c=l. 

All numerical computations within this work have been performed with Mathematica® 5.2. 



Chapter 2 

Basic Concepts of Potential Models 



This chapter presents basic concepts of quarkonium potential models, as well as a brief survey 
of modifications to the basic model introduced in recent years. 



2.1 Motivation 

Within relativistic quantum field theory the appropriate framework for the description of bound 
states is the Bethe-Salpeter formalism. In QCD with heavy charm and bottom quarks, the 
characteristic scale Aqcd ^0.2 GeV is small compared to the quark masses, rric^ 1.3 GeV and 
mf, ~ 4.5 GeV. A systematic expansion in powers of l/niq is possible (" non-relativistic QCD"). 
Bound state problems can be dealt with non-relativistically. Following these thoughts one treats 
the bound state problem by solving a Schrodinger equation using an appropriate potential. 



2.2 Schrodinger Equation 

The starting point for calculating the wavefunctions and eigenvalues of quarkonium states is the 
time-independent Schrodinger equation with a central potential and the Hamiltonian 

where we have separated the relative motion and the motion of the center of mass with 

Hem = , H^el = 2^ + ^o(rO , 

^-Pl+V2, P- rni+mz ' (22) 

7~, _ miri+m2r2 ^ _ _ ^ ^ 

mi+m2 ' ' ' 1 ' 2 ; 

M = mi + m2, mred " 



mi+m2 ' 



In the center-of-mass frame and substituting 

p ^ -iV, (2.3) 
the resulting coordinate space Schrodinger equation is: 

A 



Voir) 



ip{f) = Eip{f) . (2.4) 



2mred 

For a radially symmetric potential the wave functions written in spherical coordinates are 

i;{r,e,(j)) = Rki{r)YUd,4>) ■ (2-5) 



7 



8 



Chapter 2. Basic Concepts of Potential Models 



The radial wave functions satisfy the equation 
1 



52 2 d\ l{l + l) . , 
+ -— 1 +^ - + Vo{r) 



2wTod ydr"^ r dr J 2?7ij.ed^2 
We introduce as usual the reduced radial wavefunction 

Ukiir) = rRki{r) , 

and arrive at 



Rklir) 



EkiRklir) . 



+ + Voir) 



with the normalization condition 

d3r|V'(f)|2 = 



Ukiir) = EkiUki{r) , 



(2.6) 



(2.7) 



dndrr^{Rki{r)Y\YUO,' 



drr\Rki{r)y 



2.3 The Potential Model 



dr {ukiir)y 



1 . 



(2S 



Various quarkonium potential models have been used over the years, some of which are purely 
phenomenological e.g. |Mar80j . others use perturbative QCD as a guide for the short range 
part of the potential and a phenomenological long range part to account for confinement e.g. 
[EGK+78|[EiH79] . In the region tested by experiments most of these models coincide up to an 
adjustable constant in the energy. In this section we concentrate on the so called Coulomb-plus- 
Linear potentials, as these are the most commonly used ones. 

2.3.1 The Coulomb-plus-Linear Potential 

As the strong interaction is described by an asymptotically free theory, one expects that the 
short-distance structure of quarkonium is adequately described by perturbative QCD with a 
small "running" coupling constant. The leading term of the quark-antiquark potential - arising 
from a perturbative QCD calculation (see Appendix [B]) - is essentially Coulomb-like 



Vr 



Ojport 



( ) - 



(2.10) 



The color charges of the quark and antiquark are subject to confinement. Lattice gauge theories 
observe, in the limit of static quarks, a linear behaviour of the potential at large separation r 
(e.g. [B+OOb] ). Consequently, 



Vconiir) 



ar 



(2.11) 



is expected to be a reasonable choice for the long range part of the potential. Hence the potential 
we use in Eq. (j2.8p to obtain eigenvalues and wavefunctions for the qq bound states is given by 
the Coulomb-plus-Linear potential 



Voir) = Vb,pert(r) + Vconiir) + C 



4as 
3r 



+ ar + C. 



(2.12) 



The Coulomb-plus-linear potential (Fig. 12. ip . sometimes referred to as Cornell potential or 
funnel potential, has first been proposed by the Cornell group around Eichten [E+75[ IECtK+78 
IEGK"*"80] to reproduce the charmonium spectrum. 



2.3. The Potential Model 
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V(r) [GeV] 
2 




r[fm] 



Figure 2.1: The Coulomb-plus-Linear potential with parameters = 0.3, a = 1.0 GeV/fm. 

In contrast to many models of this type we restrict ourselves to C = in the present work. 
While there is no principle argument against a non-zero constant C, we wish to keep the number 
of parameters at a minimum and determine the leading order quarkonium masses as 



2mq + Eki , 



(2.13) 



where iJiq is the quark mass. 



2.3.2 Relativistic Corrections 

The Coulomb-plus-lincar potential itself is spin-indcpcndcnt and therefore not sufficient to re- 
produce the full structure of the quarkonium spectra. One has to introduce spin-dependent 
corrections to this potential to be able to describe the quarkonium spectra in detail. 

In the calculation of quarkonium masses, corrections to the potential Vq used in the 
Schrodinger equation. 



Vo{r)+5Vif) 



(2.14) 



which arise either from the perturbative QCD part or as corrections to the confining part of the 
potential, are usually treated perturbatively: 



Mik^^+Hj) = Mo,ki + SMkjis , 

SMkjis = J d\ V'^(r ) SV{r) ^{f) = {SV{f)) . 



(2.15) 



Level splitting are in general due to three interaction types: spin-spin, spin-orbit and tensor 
interaction. The spin-dependent potentials are expressed in terms of r, p, si and S2- One 
introduces the total spin 



Sl + S2 



(2.16) 
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the angular momentum 



and the tensor operator 



S 



12 



12 



L = fx p , 



jsi ■f){s2-f) 1 ^ 

5 o(^l ■ ^2) 



(2.17) 



(2.18) 



Experiments provide us with information on the total angular momentum J, the angular mo- 
mentum L and the total spin S of the quarkonium states. Expectation values of (si • 5*2), {L- S), 
{S12) are then expressed in terms of the eigenvalues j, I, S of total angular momentum J, angular 
momentum L and total spin S. 



Spin-spin coefficient 

We have 



S{S + 1) = (Sl + 82^ = Sl{si + 1) + S2{S2 + 1) + 2SI • S2 , 



and therefore 



Sl • S2 



S^-Slisi + l)-S2is2 + l] 



In the case of si = S2 = 1/2 this implies 



I for spin singlet S = 



Sl ■ S2 



+j for spin triplet S = 1 



(2.19) 



(2.20) 



(2.21) 



Spin-orbit coefficient 

For the calculation of {L ■ S) we employ 



J2 = j(j + i) = {L + Sf = L^ + S^ + 2{L-S) 
= l{l + l) + S{S + l) + 2{L-S), 



to find 



L-S 



1 



:i(j + !)-/(/ + !) -5(5 + 1)] . 



(2.22) 



(2.23) 



The spin-orbit term obviously vanishes for Z = or 5 = and thus only contributes for / 7^ 
and 5 = 1. Tab. 12.11 gives matrix elements of L ■ S for states with eigenvalues j, I, 5 = 1. 



Tensor coefficient 

The tensor term can be expressed via the total spin 5, using 



(5-f) 



[(Sl + S2) ■ fY = [(Si • f) + (S2 • f)f 

{si ■ f )^ + (s2 • f)'^ + 2(si • f)(s2 • r) 



1 



+ 2(si • f )(S2 • f) . 



(2.24) 



It follows that 



(si-f)(s2-r) = -{S-fy 



—r 



(2.25) 



2.4. Beyond the Basic Model 
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j 


l + l 


I 


/ - 1 


L-S 


I 


-1 


-(^ + 1) 



Table 2.1: Spin-orbit coupling for Z / and S = 1. 



With 



this leads to 



'12 



3<^ - 



It is evident that the tensor term vanishes for S = 0. Furthermore 

5i2 = for / = , 

since in this case 



A lengthy calculation gives 

Sl2 



(2Z-l)(2Z + 3) 



S^P -^L- S -3{L- Sf 



for the diagonal elements of 5*12 which are listed in Tab. 12.2 



(2.26) 

(2.27) 

(2.28) 
(2.29) 

(2.30) 



3 


l + l 


I 


/ - 1 


5*12 


21 
21+3 


2 


2(^+1) 
21-1 



Table 2.2: Non- vanishing diagonal elements of 5i2. 



2.4 Beyond the Basic Model 

The basic model is usually considered to be built of a linear confinement potential plus a po- 
tential arising from one-gluon exchange. In the attempt to describe the quarkonium spectra 
in greater detail different modifications have been applied to the model, some of them purely 
phenomenological and others of more fundamental nature. We give a brief survey over the most 
prominent ones. 
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Phenomenological as{r) 

The basic model does not include a running coupling although the short range part is calcu- 
lated using perturbative QCD. As this is a basic feature of QCD and because there is no exact 
expression for the running coupling constant as(r) in coordinate space available, a first idea is 
to model OsiQ'^) in the relevant region and transform it to coordinate space |GI85t [ZVOR95] . 

Relativistic corrections for the confinement potential 

Most of the more recent Coulomb-plus-Linear models introduce relativistic corrections to the 
confining part of their potential. These corrections are due to an assumed Lorentz-structure of 
the linear potential. Usually corrections analogous to a Lorentz-scalar [GRR82) or a mixture 
of Lorentz-scalar and Lorentz- vector exchange [EFG03[ IRROTj are considered. These terms are 
implemented as they improve spectral properties of the P-states. 

One-Gluon exchange and higher order effects 

In the basic model the one-gluon exchange with or without retardation is considered as origin of 
the short distance part of the potential. As the short distance part is due to a perturbative QCD 
calculation, higher order diagrams have been considered in some models |GRR82tlEFG031IRR07] . 



In this setting we want to study the following questions: 

• How much phenomenological input is needed to obtain a satisfactory description of the qq 
spectra? 

• How important are higher order effects for the short range part of the potential? 

• Is it possible to build a model, with minimal phenomenological input, with parameters 
close to the commonly expected values 

as{ml) « 0.3 , as{ml) 0.2 , 

mc = 1.15- 1.35 GeV, = 4.6 - 4.9 GeV , 

a w 1 GeV/fm , (2.31) 

given by QCD determinations (see E"'"04] ) or lattice simulations (e.g. B^OOaj )? 

To gain some insight into these issues, we will study Coulomb-plus-Linear models with per- 
turbatively treated corrections from various origins. To keep the phenomenological input at 
a minimum, we neither consider a phenomenological running coupling as{r) nor assume any 
Lorentz-structure, taking the confining potential to be purely static, in these models. Having 
said this, the origins of the perturbatively treated corrections in our models are either pertur- 
bative QCD (through expansion in powers of as) or relativistic corrections (through expansion 
in powers of 1/mq). 



Chapter 3 

The Breit Interaction 



This chapter motivates and discusses a potential model with perturbative corrections based on 
the Breit interaction. 



3.1 Perturbative Corrections 

The origin of the short range part of the potential is perturbative QCD which is appropriate 
for the description of high energy scattering processes. Thus we actually consider a scattering 
process, from which we deduce the potential. Potentials including relativistic corrections have 
first been deduced for electron-electron and electron-positron scattering by Gregory Breit |Bre29l 
IBre301 IBre32] . These potentials are usually referred to as Breit interaction. In analogy to 
these potentials we have rederived the Breit interaction for quark-antiquark in Appendix IB. 1^ 
which is equivalent to an one-gluon exchange potential including retardation corrections. After 
transformation into the center of mass frame, which effectively results in the substitution 



mi = m2 , {q = c, b) 
Pi = -P2 , 



(3.1) 



the Breit interaction for the quark-antiquark system (Eq. ()B.35P ) reads 



+ 



3m2 



3r 3m^ 



2as 
3m^ 



p.p if ■ p){r ■ p) 



3(si • f)(s2 ■f)-Si-S2 



(3.2) 



+ 



2as (fxp)- (si + S2) 



This is taken to be the perturbative QCD part of the potential, and thus the origin of the 
perturbatively treated corrections, used in this model. 



3.2 Potential Model 

The potential in this model is composed of a linear confining part and the Breit interaction for 
quark-antiquark (Eq. (|3.2p l: 



P ■ P _l_ {f ■ p){f ■ p) 



(3.3) 



4a. 
3m^ 



3(si • r){s2 ■f)-si-S2 



+ 



2as (fxp)- (si + S2) 



13 



14 



Chapter 3. The Breit Interaction 



As for all Coulomb-plus-Linear models, the potential used in the Schrodinger equation to find 
bound states is 

dry 

Voir) = -^ + ar, (3.4) 
Sr 

while the remaining parts of Eq. (13.30 are treated within perturbation theory. The resulting 
mass formula is given by 

^,[CpL-B](,25+i,^.) ^ 2m, + i^H + ^|^(0)P + ^Q5(5 + l)-| 

j{j + 1) -1(1 + 1) - S{S + 1) / 1\ , Su I 1 



+as 2 \ ^ 7 + " 

The last term represents the expectation value of the purely momentum dependent terms (see 
App. |D]). Though written according to option 1 in this Appendix, the expectation value for 
the momentum dependent terms in the case of the Breit interaction is unique (see App. ID.2p . 
To mark figures, tables and quantities belonging to this model they are denoted, as already 
indicated for the potential and the mass formula, with [CpL-B]. 



3.2.1 Charmonium 



To determine the model parameters we use those states for which we have reliable experimental 
data and maximum confidence in the potential model. Therefore the parameters for the model 
[CpL-B] are fixed by a fit to the masses of the states r/o J/ '4' aiid h^, resulting in the set of 
parameters 

as = 0.29 , 

rric = 1.2185 GeV, (3.6) 
a = 1.306 GeV/fm. 

One recognizes that Os is in a region we can consider as reasonable (Eq. (j2.3ip ) and that the 
charm quark mass rric is within the interval given in [E"'"04] . The value for the string tension a 
does not match with the usual value found in lattice simulations but is about 30% larger than 
expected. 

The resulting masses are compared with experimental data in Tab. 13.11 while in Fig. 13.31 
the emerging spectrum is confronted with the empirical spectrum. Taking a closer look at 
the spectrum one realizes that, for states below the DD-threshold, model [CpL-B] is able to 
reproduce the general structure, though there are a few issues concerning details of the spectrum. 
Firstly, the spin-spin splitting for the 2S'-states is too large in magnitude. This is easy to 
understand by inspecting the radial wavefunctions (Fig. 13.11 &: [312]) . as the spin-spin splitting 
in model [CpL-B] is solely depending on the wavefunction at the origin |'0fc;(O)|. This quantity 
is related via a power series ansatz to the slope of the reduced radial wavefunction Uki{r) at 
the origin. Investigating the reduced radial wavefunctions for the S-states (Fig. 13. ip . we realize 
that the slopes near the origin are almost equal for all of these states. Considering, after a quick 
glance at the radial densities in Fig. 13.41 the linear part of the potential to be dominant for the 
charmonium bound states this is not surprising as one knows that for a purely linear potential 



|V'fci5(0)| = 1^^^5(0)1, /ci,A;2GlN, 



(3.7) 



3.2. Potential Model 
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IS wavefunction 
2S wavefunction 
3S wavefunction 




r [fm] 



Figure 3.1: Charmonium S-state reduced radial wavefunctions calculated within model [CpL-B] 
{as = 0.29, a = 1.306 GeV/fm, rUc = 1.2185 GeV). 



u(r) [fm"''^] 




1.5 ^'-'2 



IP wavefunction 
2P wavefunction 
ID wavefunction 



2.5 



r [fm] 



Figure 3.2: Charmonium P- and D-state reduced radial wavefunctions calculated within model 
[CpL-B] {as = 0.29, a = 1.306 GeV/fm, = 1.2185 GeV). 
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Mass [GeV] 



4.2 - 
4 - 
3.8 - 



DD-threshold 
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3.4 - 

3.2 - 

3 - 

2.8 - 



'Xcl- 

Xco- 
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Figure 3.3: Charmonium spectra of experiment and model [CpL-B] {as = 0.29, a 
1.306 GeV/fm, nic = 1.2185 GeV); for experimental data masses and widths are shown. 




Figure 3.4: Radial densities for charmonium states of model [CpL-B] plotted together with used 
potential (a^ = 0.29, a = 1.306 GeV/fm, rric = 1.2185 GeV); base lines of radial densities have 
been shifted by according Eki- 



3.2. Potential Model 
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Experiment E+04 


TViPorv fPnl, Rl 


State 


V..' cii.j.\j.nj.cij 








Mass [MeVJ 


Widtn [iVleVJ 


Ma'^c; fMpVl 

lV±cloo IVXC V 


i OQ 




2980.4 ± 1.2 


25.5 lb 3.4 


9QSn 


1 
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U.UyiU it U.UUoz 


ouy / 


-f^ 1 
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3526.21 ± 0.25 
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3526 


1 3 P„ 
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o4io.l0 =b U.oO 


iU.Z it U.o 




l^Pi 


Aci 


0010.59 ± 0.10 


0.9d ± 0.12 


3501 


1 3 Po 




o55d.2d ± 0.11 


2.25 ± 0.15 




Z i^o 
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14 ± 7 


ouy c/ 
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U.zoo it U.Ui / 
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1 ^ Do 
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1//' 


3770 lb 2.4 


23.6 lb 2.7 


■^870 
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41 OQ 


93 P„ 










23 Pi 








4084 










4143 










4260 


3^Si 




4040 lb 10 


52 lb 10 


4347 



Table 3.1: cc-state masses from experiment and model [CpL-B] (a^ = 0.29, a = 1.306 GeV/fm, 
rric = 1.2185 GeV); masses and widths are displayed for the experiment. 



is valid. To reproduce the experimental spin-spin splitting for the 2S'-states one would need 

- 0.64 , (3.8) 
IV'is(0)| ^ ^ 

instead. As other models jEFG03t IGJ961 IFul91l IGI85] face the same problems, we assume that 
the vicinity of the DD-threshold leads to non-negligible corrections for these states. Secondly, 
one realizes that, even though the mag nitude of M{l^P2) - M{I^Pq) differs only shghtly for 
model [CpL-B] and experiment, the model does not reproduce the triplet P-state splitting 
satisfactorily. This splitting is generated by the tensor and spin-orbit parts of the Hamiltonian, 
which are given for model [CpL-B] by 

H-^^H- = (o£.5 + 5..)^(i). (3.9) 

Thus the proportion 

^[CpL-B]^p^ _ ( MiXc2)-MiXcl) \ _4 
^ ^ ^ " U(Xcl)-M(Xe0)j[CpL-B]"5' ^'-'"^ 

does not depend on parameters and will never reproduce the experimental value 

e^(l^) = (wMr^^) =0-482 ±0.005, (3.11) 

leading to the conjecture that the Coulomb-plus-Linear model [CpL-B] is not sufficient to de- 
scribe the spectrum in detail. 
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3.2.2 Bottomonium 

' """""" ""^ no 



In contrast to charmonium there are, apart from first hints for ijh measured at ALEPH H^02 
masses for bottomonium singlet states known. Thus we cannot fix the parameters of the model 
in the same manner as for charmonium. The states we use as input to determine parameters 
are T(IS'), T{2S) and C{1P), the center of gravity for the IP triplet states defined as 

C{1P) = ^(5M(x62) + 3M(xw) + M(x6o))~ 9900MeV. (3.12) 

Identifying the C{1P) with the l^Pi state of the model the resulting parameter set is 

as = 0.388 , 

TTib = 4.7645 GeV, (3.13) 
a = 1.02GeV/fm. 

The string tension a and the bottom quark mass are in the expected regions, but starting 
from the value of Os found for charmonium, one would expect a much smaller value for a^. As 
we will see this is due to the special properties of the IS states. 

The masses computed within model [CpL-B] are compared to experimental data in Tab. 13. 2^ 
while in Fig. 13.71 the resulting spectrum is displayed together with the experimental one. 







Experiment jE"'"04j 


Theory [CpL-B] 


State 


Candidate 








Mass [MeV] 


Width [MeV] 


Mass [MeV] 


1^5o 


(%) 


(9300 ± 20 ± 20) 




9283 


l^Si 


T(1S) 


9460 ± 0.26 


(53.0 ± 1.5) keV 


9460 


I'Pi 








9900 


l^Po 


XboilP) 


9859.44 ± 0.42 ± 0.31 




9850 


l^Pi 


XbiilP) 


9892.78 ±0.26 ±0.31 




9888 


1^P2 


Xb2{m 


9912.21 ±0.26 ±0.31 




9918 


2' So 








9946 


2^ Si 


T(2S) 


10023.26 ± 0.31 


(30.6 ± 2.3) keV 


10023 










10172 










10158 




T{1D) 


10161.1 ±0.6 ± 1.6 




10169 


l^Da 








10180 


2^ Pi 








10274 


2^Po 


Xbo{2P) 


10232.5 ±0.4 ±0.5 




10233 


2^ Pi 


Xbi{2P) 


10255.46 ± 0.22 ± 0.50 




10264 


2^P2 


Xb2{2P) 


10268.65 ± 0.22 ± 0.50 




10289 


3' So 








10324 


3^51 


T(3S) 


10355.2 ± 0.5 


(22.1 ± 2.7) keV 


10383 


4^50 








10626 


43 Si 


T(4S) 


10580.0 ± 3.5 


20 ± 2 ± 4 


10678 



Table 3.2: 66-state masses from experiment and model [CpL-B] (a^ = 0.388, a = 1.02 GeV/fm, 
= 4.7645 GeV); masses and widths are displayed for the experiment. 
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IS wavefunction 




2S wavefunction 




3S wavefunction 




4S wavefunction 




r [fm] 



Figure 3.5: Bottomonium 5-state reduced radial wavefunctions calculated within model [CpL-B] 
{as = 0.388, a = 1.02 GeV/fm, nib = 4.7645 GeV). 
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Figure 3.6: Bottomonium P- and I?-state reduced radial wavefunctions calculated within model 
[CpL-B] {as = 0.388, a = 1.02 GeV/fm, rrib = 4.7645 GeV). 
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Mass [GeV] 

10.6 BB-threshold 



10.4 



10.2 



10 



9.8 - 

9.6 - 

9.4 - 

9.2 - 



Experiment 

ALEPH 
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'So 



'Si 



Figure 3.7: Bottomonium spectra of experiment and model [CpL-B] (og = 0.388, a = 
1.02 GeV/fm, nib = 4.7645 GeV); ijf, mass measurement by ALEPH collaboration is included 
separately in the experimental spectrum. 




Figure 3.8: Radial densities for bottomonium states of model [CpL-B] plotted together with 
used potential (a^ = 0.388, a = 1.02 GeV/fm, nih = 4.7645 GeV); base lines of radial densities 
have been shifted by according Eki- 
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The reduced radial wavefunctions calculated within model [CpL-B] are plotted in Fig. 13.51 and 
Fig. 13.61 Examining Fig. 13.51 one realizes that the slope at the origin of the IS'-state differs 
significantely from the slopes of the rest of the S'-states, while their slopes again are all nearly 
equal. A quick glance at the radial densities (Fig. 13. 8p shows that the 15-state is lying way 
deeper down in the funnel than all the other states. Both observations emphasize the special 
role of the IS'-state for bottomonium, as it is, unlike the other states, strongly influenced by the 
Coulomb part of the potential. 

From the spectrum (Fig. 13. 7p one realizes that the model matches experimental data quite 
well, with the exception of the triplet P-state splittings. The ratios 



[CpL-B] 



4 

5 



[CpL-B]. . ^ ( M{xb2)-Mixbi) \ 
' ^ ' \M{xb,)-M{xm))\ 

' ^ ^ " U(4i)-M(x'jj.cpL-Bl"5' ^'-'"^ 



[CpL-B] 

are independent of the model parameters and do not coincide with the experimental values 
$fP(lP) = f ^;^^^)-^;^''!^ =0.58 ±0.03, 



' / Exp 
61^' 



ci,fP(2P) = =0.57 ±0.05. (3.15) 

' \M{x',,)-M{x'JJe.p 

This reenforces the presumption that the Coulomb-plus-Linear model [CpL-B] is not sufficient 
to describe quarkonium spectra below threshold in detail. On the other hand the triplet P- 
state splittings are enlarged. These deviations in the overall size of the P-state splittings are 
understandable, since the magnitude of the relativistic corrections responsible for those splittings 
are governed by the strong coupling constant and we already argued that the value one would 
expect for as is significantly smaller than the one used in this model. The agreement between the 
1^ So-state mass of model [CpL-B] and the (tentative) mass found by the ALEPH collaboration 



H"'"02 is a coincidence. 



3.3 Mass-Radius Relationship 

After studying the spectra generated by model [CpL-B], we examine whether the model leads to 
an approximately linear mass-radius relatior0. This is of interest in the study of string breaking 
[P+05l[BND+05] . 

Using the virial theorem for each given eigenstate \^/J), 

(V|r|v^) = l(V^|f-w(f)|v), (3.16) 

with T being the kinetic energy and V{f) the potential used in the Hamiltonian 

H = ^ + , (3.17) 
one finds for a radially symmetric potential: 

mij) = ^{i^\rV'{rm . (3.18) 



^Of course, in the present context, the term "radius" means average distance between quark and antiquark. 
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In the case of V{r) = —-^ + err we have 

E = {i^\H\ij) = {i;\T\i;) + {^P\V\^l^) = ^{^^\rV'\^l^) + {'^P\V\i;) 

Adding the quark masses this leads to 

M = 2mg + E = 2mg + —{r)-^(-j, (3.20) 

with the leading term proportional to (r) accompanied by a small offset due to the Coulomb 
part. 

When probing the relation of mass and radius one has to decide which mass to take into 
account as the perturbatively implemented corrections do not change the radius. One possibility 
is to consider only the bare masses, defined by 

= 2mg + Eki , (3.21) 

not including any of the perturbatively treated corrections. 

3.3.1 Charmonium 

It is instructive to have a look at the radial densities in Fig. 13.41 before investigating the mass- 
radius dependence in detail. It occurs that all the states we consider are dominated by the linear 
part of the potential and thus the linear radius- mass relation (Eq. ()3.20p ) is realized neglecting 
the small correction from the Coulomb term. 

In Tab. 13.31 bare masses and radii for states in model [CpL-B] are presented. 



State 


Mll'^ [MeV] 


(r) [fm] 


15 


3086 


0.397 


IP 


3573 


0.616 


2S 


3828 


0.752 


ID 


3944 


0.796 


2P 


4179 


0.918 


35 


4403 


1.037 



Table 3.3: Charmonium radii and bare masses in model [CpL-B] (a^ = 0.29, a = 1.306 GeV/fm, 
nic = 1.2185 GeV). 



Matching a straight line to the data in Tab. 13.31 vields for the parameters of the best fit straight 
line 

M = a + b{r) , 
a = 2.285 GeV, (3.22) 
b = 2.060 GeV/fm. 
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The resulting straight Une is displayed together with input data in Fig. 13.91 As expected the 
fitted straight line reflects the dependence of mass and radius almost perfectly, even the virial 
theorem approximation for the slope of the straight line 



dM 



d{r) 2 

holds with high accuracy. 



-a , (3.23) 



3.3.2 Bottomonium 

As for charmonium it is instructive to inspect the radial densities (Fig. 13. 8p before studying the 
mass-radius relation in detail. One realizes that for bottomonium the IS-state is dominated by 
the Coulomb-part of the potential, while the remaining states seem to be mainly dominated by 
the linear part. 

In Tab. 13.41 the radii and bare masses calculated within model [CpL-B] are presented. 



State 


M|'f " [MeV] 


(r) [fm] 


15 


9421 


0.186 


IP 


9921 


0.364 


2S 


10028 


0.448 


ID 


10185 


0.502 


2P 


10301 


0.583 


35 


10397 


0.657 


45 


10696 


0.837 



Table 3.4: Bottomonium radii and bare masses in model [CpL-B] (a^ = 0.388, a = 1.02GeV/fm, 
mft = 4.7645 GeV). 



In Fig. I3.10l the bare masses are plotted over the radii of states together with the best fit straight 
line, which parameters are given by 

M = a + b{r) , 
a = 9.336 GeV, (3.24) 
b = 1.515 GeV/fm, 

where the 15-state is not included into the fit. Though the linear behaviour for bottomonium 
is not as perfect as in the case of charmonium, the deviations of the considered states from the 
straight line are relatively small. Further one finds that the virial theorem approximation for 
the slope of the straight line 

dM 3 , , 



is again approximately fulfilled. 



24 



Chapter 3. The Breit Interaction 



Mass [GeV] 

4.75 
4.5 
4.25 
4 

3.75 
3.5 
3.25 

3_ 



♦ 


S -States 


A 


P-states 


★ 


D- states 




Fit 



0.4 



0.6 



0.8 



1.2 



<r> [fm] 



Figure 3.9: M^^f^^ versus (r) for cc-states in model [Cpl-B] (a^ = 0.29, a = 1.306 GeV/fm, 
nic = 1.2185 GeV). 
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Figure 3.10: MI^''' versus (r) for 66-states in model [Cpl-B] (a, = 0.388, a = 1.02 GeV/fm, 
= 4.7645 GeV); IS-state is not included in straight line fit. 



3.4. Summary 



25 



3.4 Summary 

In the study of model [CpL-B] we have seen that the spectra arising from the model already bear 
resemblance to the corresponding experimental spectra, though we encountered some problems. 

It occurs that the string tension for charmonium in model [CpL-B] is about 30% larger than 
the empirical value. Furthermore we find that the values for the strong coupling constant for 
bottomonium and charmonium behaves in an unusual way, giving a larger value for bottomonium 
than for charmonium. As a consequence the spin-dependent splittings in bottomonium are, with 
respect to the experimental spectrum, generally enlarged. 

For charmonium it is not possible to reproduce the spin splitting for the 2S'-states. The near 
£)D-threshold may have a non-negligible influence on these states. The most severe problem 
of model [CpL-B] is that the experiment provides a different structure for the triplet P-state 
splitting. The model yields the built-in value ^\nP) = 4/5 independently of parameters 

or radial excitations, while the experiment significantly differs from this value both for bottomo- 
nium and charmonium. Based on this observation we conclude that the discussed model may 
be a first approximation for a quarkonium potential model, but has to be improved in order to 
reproduce details of the quarkonium spectra. 

In addition we investigated the mass-radius dependence for states of charmonium and bot- 
tomonium. We confirmed a leading linear relation between mass and radius for charmonium and 
found that, with the exception of the IS-state, the linear relation is also a good approximation 
for bottomonium. 



Chapter 4 



One-Gluon Exchange 
without Retardation 



This chapter motivates and discusses a potential model with perturbative corrections based on 
One-Gluon exchange without retardation. 

4.1 Perturbative Corrections 



In Chapter [3] we have investigated a Coulomb-plus-Linear model with perturbative corrections 
based on Breit interaction, model [CpL-B]. In the calculation of the Breit interaction (see App. 
IB.ip we explicitly used the on-shell condition for the energy transfer, 



E'a — Ea — Eb — E't 



B > 



(4.1) 



in the expansion of the propagator (Eq. (|B.30|) ). leading to retardation. While suitable for the 
case of scattering, using this condition in the case of bound states is not valid. For a bound 
state (7o = should be the more natural choice. The corresponding potential (see App. IB.2P is 
the one arising from one-gluon exchange without retardation corrections. Transformed into the 
center of mass frame, and using Eq. dM]), the potential Eq. (IRM]) reads 



y{2)( 



6r drrig 



Aas p- p , 2as {r x p) ■ {si + S2) 



Snig r 



+ 



m„ 



+ 



Sl ■ S2 



(4.2) 



4.2 Potential Model 



The quark-antiquark potential in this model is composed of a linear confining part and the 
one-gluon exchange potential without retardation (Eq. (j4.2p l: 



Sr Smi 



4as p- p , 2as {f x p) ■ (si + S2) 



3m^ r 



+ 



+ 



(4.3) 



As for all Coulomb-plus-Linear models, the potential used in the Schrodinger equation to find 
bound states is 



Voir) 



4a, 
3r 



-|- ar 



(4.4) 
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while the remaining parts of Eq. (14.30 are treated perturbatively. The resulting mass formula 
is given by 



2 



j{j + l)-l{l + l)-S{S + l) / 1 \ , Su/l 



+^2(^<^'^^*inl{v'^{n)-^d\l\V^Pif)\A . (4.5) 



3m'^ \ 

Sm^ \J ' ' r 

The last term represents the expectation value of the purely momentum dependent terms using 
option 3 in Appendix [Dj Figures, tables and quantities related to this model are denoted as 
[CpL-3]. 

4.2.1 Charmonium 

The parameters in model [CpL-3] are determined as in model [CpL-B], namely by a fit to the 
masses of states r/c, J j-il) and /ic, providing the parameter set 

a, = 0.302 , 

rUc = 1.272 GeV, (4.6) 
a = 1.21 GeV/fm . 







Experiment [E+04] 


Theory [CpL-3] 


State 


Candidate 








Mass [MeV] 


Width [MeV] 


Mass [MeV] 


I'So 




2980.4 ± 1.2 


25.5 ±3.4 


2980 






3096.916 ±0.011 


0.0910 ±0.0032 


3097 




he 


3526.21 ±0.25 


< 1.1 


3527 


l3Po 


XcO 


3415.16 ±0.35 


10.2 ±0.8 


3430 




Xcl 


3510.59 ±0.10 


0.96 ±0.12 


3503 


l3P2 


Xc2 


3556.26 ±0.11 


2.25 ±0.15 


3560 


2^50 


ic 


3638 ± 5 


14 ±7 


3674 


23^1 


^' 


3686.093 ± 0.034 


0.283 ±0.017 


3765 










3891 


l^Di 




3770 ± 2.4 


23.6 ±2.7 


3855 










3884 


13^3 








3911 


2iPi 








4077 


23Po 








3985 


23 Pi 








4054 


23P2 








4109 










4207 


33^1 




4040 ± 10 


52± 10 


4291 



Table 4.1: cc-state masses from experiment and model [CpL-3] (a^ = 0.302, a = 1.21 GeV/fm, 
rric = 1.272 GeV); masses and widths are displayed for the experiment. 
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The resulting parameter values for quark mass nic and strong coupling constant as are similar 
to those in model [CpL-B]. The string tension a decreased by about 10% with respect to model 
[CpL-B] and is somewhat closer to the value suggested by lattice simulations. 

The results of the model calculations are presented in comparison with experimental data in 
Tab. 14.11 In Fig. 14. li the spectrum is displayed together with the cc-spectrum from experiments. 
As model [CpL-3] differs form model [CpL-B] only in the purely momentum dependent term, 
model [CpL-3] faces the same problems for the spin-dependent splittings as model [CpL-B], 
namely the built-in ratio 

^[CpL-3] ^ ( M{xc2)-M{x.^) \ ^4 
\M{xci)- M{xco)J 

[CpL-3] ^ 

which is not consistent with experimental data (Eq. (|3.11|) ). The decrease of the string tension 
a compared to model [CpL-B] leads to a slight reduction in size of the triplet P-state splitting. 



4.2.2 Bottomonium 

The parameters of model [CpL-3] 



as = 0.326 , 

mb = 4.7215 GeV, (4.8) 
a = 1.12GeV/fm, 







Experiment jE"'"04j 


Theory [CpL-3] 


State 


Candidate 








Mass [MeV] 


Width [MeV] 


Mass [MeV] 


l^^o 


(%) 


(9300 ± 20 ± 20) 




9351 


l^^i 


T{1S) 


9460 ± 0.26 


(53.0 ± 1.5) keV 


9460 










9901 


l3Po 




9859.44 ± 0.42 ± 0.31 




9861 




XbiilP) 


9892.78 ±0.26 ±0.31 




9890 




Xb2{lP) 


9912.21 ±0.26 ±0.31 




9915 


2' So 








9967 


2^ Si 


T{2S) 


10023.26 ± 0.31 


(30.6 ± 2.3) keV 


10023 










10170 










10158 




T{1D) 


10161.1 ±0.6 ± 1.6 




10168 










10177 


2^ Pi 








10284 


2^Po 


Xm{2P) 


10232.5 ±0.4 ±0.5 




10250 


2^ Pi 


Xbi{2P) 


10255.46 ±0.22 ±0.50 




10276 


2^P2 


Xb2{2P) 


10268.65 ±0.22 ±0.50 




10296 


^'So 








10350 


3^ Si 


T(35) 


10355.2 ± 0.5 


(22.1 ± 2.7) keV 


10395 










10663 


43^1 


T(4S) 


10580.0 ± 3.5 


20 ± 2 ± 4 


10703 



Table 4.2: 66-state masses from experiment and model [CpL-3] [as = 0.326, a = 1.12 GeV/fm, 
rrih = 4.7215 GeV); masses and widths are displayed for the experiment. 
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DD-threshold 
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he 
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Figure 4.1: Charmonium spectra for experiment and model [CpL-3] (a^ = 0.302, a 
1.21 GeV/fm, mc = 1.272 GeV); for experimental data masses and widths are shown. 
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Figure 4.2: Bottomonium spectra from experiment and model [CpL-3] (a^ = 0.326, a = 
1.12 GeV/fm, nib = 4.7215 GeV); % mass measurement by ALEPH collaboration is included 
separately in the experimental spectrum. 
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are determined in the same manner as in model [CpL-B], by a fit to the masses of states T(15), 
T(25) and C(1P) (Eq. UjSJTh ). identifying C{IP) with the l^Pi state of the model. The 
strong coupling constant in model [CpL-3] is now significantely smaller than in model [CpL-B] 
but still definitely too large with respect to the strong coupling constant of model [CpL-3] for 
charmonium. Furthermore the reduction in comes at the cost of a slightly increased string 
tension, while the quark mass stays about the same. 

The resulting masses are compared to experimental data in Tab. 14. 2[ In Fig. 14.21 the 
emerging spectrum is compared to the experimental one. The spin-dependent parts of [CpL-3] 
and [CpL-B] are identical. The problem with the triplet P-state splitting prevails, resulting once 
more in the ratios 



[CpL-3]. . ^ f M{xb2) - Mixbi) \ ^4 



4 

V^(4i)-M(x',o)y[CpL-3]^5' ^^-^^ 

which disagree with the experimental values (Eq. (|3.15|) ). Due to the smaller Og the total size 
of the triplet P-state splitting is now close to the empirical one. 

4.3 Summary 

The study of model [CpL-3] is mainly inspired by the fact that the on-shell condition used in 
the computation of the Breit interaction is not valid for the investigated system. 

As model [CpL-3] differs only in the purely momentum dependent terms from model [CpL-B], 
it provides similar spectra for the spin-dependent splittings as model [CpL-B]. Problems remain 
for the triplet P-states and to a lesser degree for the 25-state spin-spin splitting in charmonium. 

The problem with the 2S-states for charmonium seems to be due to the near P'D-threshold, 
which we do not account for in any of the models. 

The momentum dependent terms in [CpL-3] are not uniquely determined. The spectra 
depend on the ordering of gradient operators. 

Model [CpL-3] yields, compared to model [CpL-B], a slightly decreased string tension for 
charmonium, and a significantly decreased but still definitely too large strong coupling constant 
Us for bottomonium. 



Chapter 5 



One-Gluon Exchange plus Induced 
Interaction 



The first section of this chapter is dedicated to the motivation and computation of the induced 
interaction, which is used to improve the triplet P-state structure of the Coulomb-plus-Linear 
model. The remaining sections present the improved model and discuss its results. 

5.1 Induced Interaction 

In Chapters [3] & m we observed that the Coulomb-plus-Linear models, with perturbative correc- 
tions arising from perturbative processes to Order 0{as), are not sufficient to describe the triplet 
P-state structure correctly. We now take a closer look at the Schrodinger equation, which we 
used to determine eigensolutions and eigenvalues. The potential used in the Schrodinger equa- 
tion is composed of a phenomenological confinement part and a short range part derived from 
perturbative QCD. The corresponding Feynman diagrams are shown in Fig. 15.11 where the 
one-gluon exchange diagram is introduced as a symbol to represent the perturbative QCD part. 



Figure 5.1: Perturbative QCD short range part and phenomenological confining part composing 
the potential used in the Schrodinger equation. 

Solving the Schrodinger equation is equivalent to solving the Lippmann-Schwinger integral 
equation [Lan90] 




IV') 



\^) + GV\(t>) + GVGV\<p) + GVGVGV\<p) + ... . 



(5.1) 



With the definition of the T-matrix 



TE{k',k) 



(5.2) 
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the Lippman-Schwinger equation for the T matrix is 

{c|>^'\V\^l^^) = {^^,\V\<t>^) + {cl)^,\VGV\i;j;) , (5.3) 

or using Eq. (jS.ip and rewritten as an operator equation: 

T = V + VGV + VGVGV + ... . (5.4) 

The Lippmann-Schwinger is commonly used for scattering problems. Bound states are identified 
as Poles in the T matrix Te at energies -E < 0. 




Figure 5.2: Diagrams corresponding to the operator equation of the T- matrix. 

One realizes that the Schrodinger equation iterates the interaction to all orders and that finding 
an eigensolution corresponds to the formation of a bound state. The formation of a bound state 




Figure 5.3: s-channel diagram corresponding to the formation of a bound state. 

is equivalent to the s-channel diagram in Fig. 15.31 Crossing symmetry implies that, for any 
s-channel process forming a bound state of a particle-antiparticle system as in Fig. 15. 3^ there 
is also a t-channel process as shown in Fig. 15.41 This corresponds to the exchange of a bound 
state between quark and antiquark. 

In other words, by solving the Schrodinger equation we find bound states which correspond 
to s-channel processes. These s-channel processes induce t-channel interactions which have to 
be taken into account. While the s-channel diagram involves a sum over ladders (Fig. 15. 2p . the 




Figure 5.4: t-channel diagram corresponding to the exchange of a bound state. 
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t-channel diagram represents the summation of bubbles (Fig. I5.5p . Solving a Schrodinger or 
Bethe-Salpeter equation yields an approximation based on ladder diagrams, whereas taking into 
account the induced interaction generates a solution based on a more complete mix of bubble 
and ladder diagrams. 




Figure 5.5: Summation of bubbles corresponding to the exchange of a bound state. 



The Lorentz-structure of these additional interactions are given by the properties of the 
corresponding bound states. In principle one would have to take into account all bound states, 
but for our model we consider only the IS'-states. These states should have the strongest 
influence on the spectrum, as they have the lowest mass and the effective range of the interaction 
decreases with increasing mass of the exchanged bound state. The IS-states for quarkonium 
are a pseudoscalar meson and a vector meson, hence we have to consider a pseudoscalar and a 
vector particle exchange interaction. 



5.1.1 Pseudoscalar Meson Exchange 

The amplitude in Feynman gauge corresponding to the massive pseudoscalar exchange diagram 
(Fig. 15. 6p is, utilizing that initial and final quark-antiquark states are color neutral, given by 

Mfi = [u{p'a, s'a) (ifi-pTs) u{PA, Sa)] 2 ^ A/r2 K^'-B' {}gplb) v{p'b, s'b)] 

q ivip 

- [u{p'A,s'Ah5y'{PA,SA)v{pB,SBh5v{p'B,s'B)] ■ (5.5) 



g2-M2 

Applying the normalization conventions as in Appendix lA.H the Dirac spinors read 



u{pa,sa) 



' EA + ruA I ^ XsA 
2mA \ Eif^Xs^ 



,c 



-, \ Eb + rUB ( c\ a-pB cf 

v{pb,sb) = Y 2mB V ^'^E^^ -Xs\ 



with energies Ex = \ Px ~^ "^x quark masses mx ■ Employing Eq. ()5.6p , the pseudoscalar 
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qi{pA,SA) 




QiiPB,SB) 



Figure 5.6: Quark-Antiquark pair exchanging a pseudoscalar meson. 



Dirac currents for particles and antiparticles are given by 

t a.p- \ / 1 



uip'A,s'A)75uipA,SA) 



X 



t 



a -PA 
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^■Pa 
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Xs 
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x: 



X. 



t ( f^-PB 

SB 



<^-Pb 



Eb + mB E'^ + uibJ 



(5.7) 



(5.8) 



with Pauhspinors Xs- Putting things together, remembering the definition of spin for particles 
(Eq. (IB.16P ) and antiparticles (Eq. ()B.22j) ). yields for the amplitude 



M 



fi'- 



9l 



cr -PA 
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9p (si ■ q) {h ■ q) 



where we have used 



mim2 
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1 90=0 



g2 - M2 q^ + Ml 



+ O (m^^) , (5.9) 



(5.10) 



The leading term in the amplitude (Eq. (j5.9p ) is of order 0(?77 ^), hence the relation between 
potential and amplitude 



1 



(2vr)^ 



l + 0[ — ) \M, 

777^ 



(5.11) 



does not provide any further terms, as we are only expanding up to order 0(777. ^). The Fourier 
transformation back into coordinate space 



V{r-pi,p2) = j d?qe-''l-''V{q 



(5.12) 



of the terms in Eq. ()5.9p are given by 



(27r)' 



(Tqe 



-iq-r 



g-Mpr 



q^+Ml Airr ' 



(5.13) 



d q qiqj e' 



-iqr 



q'^ + Ml 



WH --Wi 



5.1. Induced Interaction 



35 



with 
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(5.14) 



Collecting the pieces we eventually arrive at the potential for a massive pseudoscalar exchange 
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A7rmim2 
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Ml H ^ + 
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5.1.2 Vector Meson Exchange 

Employing Feynman rules one derives the amplitude in Feynman gauge for the massive vector 
exchange diagram (Fig. 15. 7p : 



M 



fi 



72 - M? 



,2 - M2 

Npa; SA)7''^i(PA, SA)t^(PB, sb)1i,v{p'b, s'b)] 



(5.16) 



It is obvious that the amplitude for the massive vector exchange (Eq. (|5.16|) ) is similar to 
the amplitude for the one-gluon exchange (Eq. ()B.6P ). apart from the color factor and the 
propagator. One can follow the computation of the potential for the one-gluon exchange, with 
the replacement 



3g2 



q2 - M2 ' 



(5.17) 



yielding 
1 



(2vr)^ 



+ ■ 



K(g;pi,P2) = 
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(5.18) 
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4mim2 

Using the equivalent to Eq. (|5.1U|) for the expansion of the propagator 

1 go=0 1 
g2 - M2 ~ g 2 + ' 

the momentum space potential for the massive vector exchange is given by 



(5.19) 
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(1i{p'b^s'b) 



qi{PA,SA) ^ ^ qj{PB,SB) 

Figure 5.7: Quark-Antiquark pair exchanging a vector meson. 



The Fourier transformation back into coordinate space 
of the various terms in Eq. (j5.20p is given by 



(5.21) 
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Collecting all the pieces we eventually arrive at the potential in coordinate space 
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5.1.3 IS-State Exchange Induced Interaction 

The potential for the 15'-state t-channel exchange interaction, referred to as induced interaction 
from here onward, is composed of the pseudoscalar meson exchange potential (Eq. (|5.15|) ) and 
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the vector meson exchange potential (Eq. (I5.24p ). Merging these potentials yields, in principle, 
two different couplings gp and g^, and two different masses for the exchanged particles Mp and 
My. The couplings gp, g^ are related to the asymptotic normalization constant of the considered 
bound state (see App. 12 of |EW88j ) . In the Coulomb-plus-Linear model the wavefunctions for 
the pseudoscalar state r]c and the vector state J /ip are identical as the potential Vq is spin- 
independent. Therefore the couplings gp and g^ are equal and we define, analogous to the 
strong coupling constant, the coupling constant for the induced interaction via 



9l 



9l 



(5.25) 



Based on the same argument, we use the degenerate IS" mass. In the case of charmonium this 
mass is taken to be the mass of a 15-state without spin-dependent splittings 



M 



1 



M{t]c) + 3M( J/V^) ) w 3.07 GeV , 



(5.26) 



while for bottomonium, due to the lack of information about the singlet mass, we choose 

M = 2mb, (5.27) 

as this choice is usually close to the mass of T(IS'). With these preparations, the potential of 
the induced interaction reads: 
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(5.28) 



(fxpi) ■ S2 - (f X P2) • si 



5.2 Perturbative Corrections 



The origin of the perturbative corrections in this model is on the one hand the one-gluon exchange 
potential without retardation (Eq. ()4.2p ) and on the other hand the induced interaction discussed 
in the first part of this chapter. Transformed into the center of mass frame by utilizing the 
substitution Eq. (|3.1|) the potential for the induced interaction reads 
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(5.29) 



By analysing this potential one notices that the perturbative corrections provided by the induced 
interaction differ in structure from the ones obtained by the one-gluon exchange (Eq. ()4.2p ). In 
particular, the relative strength of the tensor and spin-orbit terms and their radius dependence 
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differ from the ones in one-gluon exchange. The induced interaction might therefore be able 
to provide, together with the perturbative QCD part, a satisfactorily description of the triplet 
P-state splittings. 

We have not yet succeeded in establishing a rigorous relation between the asymptotic nor- 
malization constant and the coupling constant of the induced interaction. This is mainly due to 
the linearly rising potential. There exists no upper boundary and therefore no binding energy. 
In practice we treat as free parameter in this model. 



5.3 Potential Model 

The potential in this model is composed of a linear confining part, the one-gluon exchange 
potential (Eq. (|4.2p l and the potential for the induced interaction (Eq. (|5.29|) ): 
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As for all Coulomb-plus-Linear models the potential used in the Schrodinger equation to find 
bound states is 



4a, 
3r 



+ err , 



(5.31) 



while the remaining parts of Eq. ()5.30p are treated within perturbation theory. The resulting 
formula for the mass of quarkonium states is given by 
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The last term represents the expectation value of the purely momentum dependent terms with 
option 3 (see App. [D]). Figures, tables and quantities related to this model are marked as 
[CpLpI-3]. 



5.3.1 Charmonium 

The determination of parameters in model [CpLpI-3] is accomplished by a fit to the masses of 
states r/c, J/ip and he, and the ratio ^^^^{IP) (Eq. (|3.11|) ). with the additional constraint that 
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the parameters a^, rric and a should be as close as possible to their empirical values. This leads 
to the parameter set 

as = 0.25 , 
Oi = 1.1 , 

rUc = 1.2877 GeV, (5.33) 
a = l.llGeV/fm. 

The charm quark mass of model [CpLpI-3] is similar to the ones obtain in the previously discussed 
models. The strong coupling constant takes a small value, while the string tension in this model 
is closer to the empirical lattice value than in the model without induced interaction. 

The masses calculated in model [CpLpI-3] are compared with experimental data in Tab. 15.11 
The emerging spectrum is displayed together with the spectrum deduced from experiment in 
Fig. EEl 

The main motivation to implement the induced interaction has been to reproduce the triplet 
P-state structure in greater detail. The resulting model ratio 

$[CpLpi-3] ^,p^^( « 0.5 , (5.34) 

Vm(xc,i)-M(xc,o);[CpLpI-3] 

is now very close to the experimental value 

^>^^P(1P) = 0.482 ± 0.05 . (5.35) 







Experiment E+04j 


Theory [CpLpI-3] 


State 


Candidate 








Mass [MeV] 


Width [MeV] 


Mass [MeV] 




Vc 


2980.4 ± 1.2 


25.5 ±3.4 


2980 






3096.916 ±0.011 


0.0910 ± 0.0032 


3097 




he 


3526.21 ±0.25 


< 1.1 


3526 




XcO 


3415.16 ±0.35 


10.2 ±0.8 


3403 




Xcl 


3510.59 ±0.10 


0.96 ±0.12 


3505 




Xc2 


3556.26 ±0.11 


2.25 ±0.15 


3556 


2' So 


Vc 


3638 ± 5 


14 ±7 


3628 




^' 


3686.093 ± 0.034 


0.283 ±0.017 


3731 










3867 




ij;" 


3770 ± 2.4 


23.6 ±2.7 


3838 










3861 










3882 


2^ Pi 








4039 


2^Po 








3899 


2^ Pi 








4017 


2^P2 








4071 










4133 


3^Si 




4040 ± 10 


52± 10 


4232 



Table 5.1: cc-state masses from experiment and model [CpLpI-3] (a^ = 0.25, Oj = 1.1, a = 
1.11 GeV/fm, rric = 1.2877 GeV, M = 3.07 GeV); experimental widths are also displayed. 
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Mass [GeV] 



Xc2 



Experiment 
Model 



IPl 3p. 3d. 

Figure 5.8: Charmonium spectra for experiment and model [CpLpI-3] (q^ = 0.25, = 1.1, 
a = 1.11 GeV/fm, rric = 1.2877 GeV, M = 3.07 GeV); for experimental data masses and widths 
are shown. 

Although the magnitude of the triplet P-state splitting is increased by about 10%, the consid- 
eration of the induced interacton within model [CpLpI-3] leads to an enormous improvement 
compared to model [CpL-3], with its built-in ratio 

^ 5 

5.3.2 Bottomonium 

The input used to determine the parameters of model [CpLpI-3] is given by the states T(IS'), 
T{2S) and C(1P), and the ratio ^>^^^(1P), with the additional constraint that the parameters 
Us, mi, and a should be as close as possible to their empirical values. In the analysis of model 
[CpLpI-3] we found that we are able to improve the structure the triplet P-state splitting, but 
this improvement comes at the cost of parameters which are even further apart the empirical 
values as the ones in the model without induced interaction, or the incapability to match the 
states T(IS'), T(2S) and C{1P). It occurs that this is again a problem due to the special role 
of the 15-state. 

5.4 Summary 

In this chapter we exploited crossing symmetry to obtain an induced interaction which is equiv- 
alent to the t-channel exchange of bound states. By considering only the lowest lying states in 
the spectra, we computed a potential for the t-channel exchange of IS-states which arises from 
pseudoscalar and vector meson exchange processes. 

On the basis of these concepts we introduced the induced interaction to the Coulomb-plus- 
Linear model [CpL-3] , leading to additional perturbative corrections which allowed us to repro- 
duce the triplet P-state structure. 
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The model works very well for charmonium, yielding a value for the string tension close to 
the empirical one for the first time. For bottomonium we face problems due to the distinct 
influence of the Coulomb-part of the potential on the 15-states as compared to the rest of the 
bound states. The situation for the strong coupling constant does not improve up to this point, 
still taking a larger value for bottomonium than for charmonium. 



Chapter 6 



Fourth-order Gluonic Potential 



In this chapter we adopt a gluon exchange potential derived by Gupta et al., which includes 
two-gluon exchange diagrams. The potential is used, with slight modifications, as short range 
part in a Coulomb-plus-linear model. After the dicussion and motivation of the model in the 
first two sections, we examine the effect of additional corrections found by Titard and Yndurain. 



6.1 Perturbative Corrections 



The Coulomb-plus-Linear models we have discussed until now have employed an one-gluon 
exchange process to construct the short range part of the potential. The natural next step for 
a further improvement of the model is to consider not only an one-gluon exchange process but 
processes of higher order in perturbative QCD. Such a potential has already been calculated 
by Gupta et al. in the early 80's |GR81t [GR82a] . They derived a fourth-order gluonic quark- 
antiquark potential, which has later been slightly corrected by Titard and Yndurain |TY94] . 
Gupta et al. used this potential as short range part of a quarkonium potential model |GRR82) . 
The perturbative part of this potential model is, with a slight modification, adopted by us and 
taken as the short range part of the model potential. 

The modification we apply to this potential is given by a reduction of the spin-independent 
contact term and its radiative correction by a factor of two. We think that this adjustment is 
necessary as we already know the strength of the leading term of the contact structure from 
one-gluon exchange potentials. The consideration of higher order processes in perturbative QCD 
should not change the strength of the leading term. We can compare the potential in |GRR82j 
with Eq. ()4.2p and realize that the strength has to be reduced by a factor of two. The reason 
for this discrepancy seems to be the choice of kinematics. As we point out in Appendix IB. 31 
the choice of kinematics is important, meaning that different choices may lead to differences in 
the momentum dependent parts of these potentials. The kinematics choice we use to obtain 
a potential is dictated by invariance under time reversal, while the kinematics employed in 
|GR81t ITY94) lead to a potential which is not invariant under the time reversal operation (see 
App. IB.3p . In the center of mass frame the additional non-time reversal invariant terms of these 
potentials can be transformed into an additional contact term, by exploiting certain identities 
in the momentum space representation of the potential. Hence we modify the strength of the 
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contact term in (GRR82| to half of its value and obtain 
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This perturbative QCD part of the potential, to order a^, is then treated as a perturbative 
correction to the confining potential. 



6.2 Potential Model 

The potential in this model is composed of the potential given in Eq. (|6.1|) and a linear confining 
part: 



y [CpL*-2] ^) ^ ^(4) ^) + , 



(6.2) 



When trying to implement the radiative corrections to the leading Coulomb term into the 
potential used in the Schrodinger equation, this would lead to problems. One realizes that 
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(6.3) 



due to the logarithmic correction. This would imply that all wavefunctions should vanish at the 
origin 



\Mo)\ = 



(6.4) 



leading to an almost vanishing splitting for the ^Si and ^S'o-states. A meaningful procedure is 
to start from solving the Schrodinger equation with the potential 



4a, 



3r 



(6.5) 



while the remaining parts of Eq. (|6.2|) are treated as a perturbation. The resulting mass formula 
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is given by 
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The last line represents the expectation value of the purely momentum dependent term using 
option 2 (see App. |D]), as this scheme provides parameters closest to the empirical values. 

One recognizes that Eq. ()6.6p includes the number of light flavours and a renormalization 
scale. These quantities will not be treated as free parameters. Moreover the renormalization 
scale is given for a scheme by Gupta and Radford [GR82b] . which differs from the common 
MS-scheme. Gupta and Radford give the relation 
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Mms = /^GR exp 



2(33 - 2n/) 



(6.7) 



which connects their subtraction scheme with the renormalization scale of the MS scheme. We 
have checked and confirmed this relation by comparing |TY94| and |GRR82j . The number of 
light flavours is set to 



Uf 
Uf 



3 for Charmonium, 

4 for Bottomonium. 



(6., 



Usually the first guess for the renormalization scale in the modified minimal subtraction scheme 
MS would be = ruq. Titard and Yndurain maintain in [TY94J that the natural scale is 
not = mq, but determined by the average momentum of the quarks. Typically the average 
momentum for the ground state in our model is 

^p2^i/2 ^ OAnic for Charmonium , 

^^•2^1/2 0.25mfe for Bottomonium . (6.9) 
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We take 

I^MS ^ 2^'^ Charmonium , 

'"ms ~ —^b for Bottomonium . (6.10) 
3 

The analysis in |TY94j showed that this should be a good choice for bottomonium, while it might 
be unreliable for charmonium. The coupling constant evaluated at this scale is ag {ml/A) > 0.5, 
and therefore the calculation of the short range part via perturbative QCD is not justified. 
Nevertheless we consider this value for the scale as we do not know how to fix the scale more 
reliably. The model calculations are performed in the renormalization scheme by Gupta and 
Radford |GR82b| . The relation Eq. (|6.7|) results for the considered cases of nj = 3 or nj = 4 in 

/^MS~^- (6-11) 



This implies 

Charmonium: fJ-GR = it^-c , 

2 

Bottomonium: /iQR = • (6.12) 
Figures, tables, and quantities related to this model are marked as [CpL*-2]. 

6.2.1 Charmonium 

We tried to fix the parameters of model [CpL*-2] analogously to models [CpL-B] and [CpL- 
3], by matching the masses of states 7]c, J/ip^ and /ic, with the additional constraint that the 
parameters should be as close as possible to their empirical values. Problems were encountered 
in reproducing the 15 spin-spin splitting, which is, for parameters within reasonable intervals, 
significantly to small in this model. The parameter set 

as = 0.27 , 

rric = 1.326 GeV, (6.13) 
a = 1.2GeV/fm, 

is constrained by reproduction of /ic, Xi ^-iid the center of gravity IS'-state mass 

C{IS) = \{M(r^,) + 3M(J/V^)) , (6.14) 

with the spin-spin splitting as close as possible to the experimental one. 

The resulting masses are shown together with experimental values for the masses and widths 
in Tab. 16.11 In Fig l6. li the arising spectrum is opposed to the experimental one. One recognizes 
that this model yields only about 60% of the experimental 1^5i — l^^o splitting. Examining Eq. 
(16.61) one finds that there is no built-in value 



/ M(Xg,2)-M(Xg.l) \ 4 



VM(xg,i)-M(xg,o)y 5' ^^-^^^ 

as in the models [CpL-B] and [CpL-3] anymore. In model [CpL*-2] this ratio is now parameter 
dependent. With the parameter set Eq. (j6.13p one finds 

$[CpL*-2](ip) ^ f M}^^^^\Z^S^\ ^ 0.74 . (6.16) 

This illustrates, that while this proportion differs only slightly from its value for the basic 
Coulomb-plus-Linear model, it changes at least in the right direction. 
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4098 










4224 


3351 
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Table 6.1: cc-state masses from experiment and model [CpL*-2] (a^ = 0.27, a = 1.2 GeV/fm, 
rric = 1.326 GeV, /Uqr = nj = 3); masses and widths are displayed for the experiment. 



6.2.2 Bottomonium 



The parameters of model [CpL*-2] for bottomonium 

Us = 0.27 , 

nib = 4.7877 GeV, (6.17) 
a = 1.05GeV/fm, 

are determined by a fit to the masses of states T(IS'), T(25) and C(1P) (Eq. (I3.12|) ). by 
identifying C{1P) with the l^Pi state of model [CpL*-2]. Remarkably the model yields for the 
first time a coupling strength which is smaller than = 0.3, while the string tension and the 
bottom quark mass coincide with empirical values. 

The results for Eq. (16. 6p are presented together with experimental data for masses and widths 
in Tab. 16.21 while the according spectra are shown in Fig. 16.21 Examining the spectrum one 
immediately notes that the triplet P-state splitting is to small by about 25%. The proportions 



\M{xbl) - M{xbo) J [CpL*-2] 



show the same slight decrease compared to Eq. (|6.15p . which one has already observed for 
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Experiment 
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Figure 6.1: Charmonium spectra for experiment and model [CpL*-2] (a^ = 0.27, a = 
1.2GeV/fm, rric = 1.326 GeV, //qr = nT-c n-f = 3); for experimental data masses and widths are 
shown. 
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Figure 6.2: Bottomonium spectra from experiment and model [CpL*-2] (a^ = 0.27, a = 
1.05 GeV/fm, mb = 4.7877 GeV, /^gr = 2/3mb, Uf = 4); r]h mass measurement by ALEPH 
collaboration is included separately in the experimental spectrum. 
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10268.65 ± 0.22 ± 0.50 




1 0989 


3^5*0 








10371 


33^1 


T(35) 


10355.2 ±0.5 


(22.1 ±2.7) keV 


10391 










10669 


43^1 


T(45) 


10580.0 ±3.5 


20 ± 2 ± 4 


10688 



Table 6.2: 66-state masses from experiment and model [CpL*-2] [ug = 0.27, a = 1.05 GeV/fm, 
TTifc = 4.7877 GeV, /xgr = 2/3mfc, = 4); masses and widths are displayed for the experiment. 



charmonium. 

6.3 Additional Spin-Independent Terms of Order 0{m~^) 

In |TY94j Titard and Yndurain show additional terms to the potential by Gupta et al. |GRR82) . 
They maintain that the order 0{m^'^) for the spin-independent terms is not complete, plus they 
show a recipe how they calculated the part involving a fictive gluon mass. Having implemented 
perturbative corrections according to the potential computed by Gupta et al. we have to check 
the effect of these additional corrections. 

The additional terms shown in |TY94j are 

,Y/^N / 133. \q\ 28. A 128 



with A the fictive gluon mass. Titard and Yndurain define the spin-spin part in a way different 
from Gupta et al., leading to different terms in spin-independent 0(m~^). We first have to 
relate these terms to the ones we use. We subtract the spin-independent part of 



''5 + ln2-71nl2l) (2S'-3) , (6.21) 



3m2 y9 niq 
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\i (r )) [ivievj 
wiLii ii/Cj . tu.^^n 


\i (r )) [ivievj 

WlLll H/tj . t|U.i^Oj) 




1 7 A 


1 1 




li.O 


— i.O 


35 


9.5 


-2.4 


IP 


1.9 


1.9 


2P 


1.8 


1.8 


ID 


0.5 


0.5 



Table 6.3: Charmonium expectation values for the additional terms presented in [TY94j cal- 
culated within model [CpL*-2] (a^ = 0.27, a = 1.2 GeV/fm, = 1.326 GeV, hqr = rUc, 
Uf = 3). 



State 


(y(f)) [MeV] 


(y(r)) [MeV] 




with Eq. (|6.24|) 


with Eq. (|6.25|1 


15 


14.7 


3.5 


25 


7.6 


1.3 


35 


6.0 


0.8 


45 


5.2 


0.5 


IP 


0.6 


0.6 


2P 


0.6 


0.6 


ID 


0.2 


0.2 



Table 6.4: Bottomonium expectation values for the additional terms presented in [TY94| calcu- 
lated within model [CpL*-2] (a, = 0.27, a = 1.05 GeV/fm, ruh = 4.7877 GeV, /icR = 2/3mh, 
Uf = 4). 



from Eq. (I6.20p to obtain the additional terms in our definition: 

4a1 



3m2 



II In a 

18 rrin 



28, A 36 1 

— In \ h-ln2 

9 m„ 5 2 



Transforming to r-space this gives 



y(f) 



4a. 
3m2 



77 ^2 /^ln(mgr) +7E 
727r 



\ /^28, A 36 1, \ 



(6.22) 



(6.23) 



It is obvious that the fictive gluon mass A poses a problem for the calculation of an expectation 
value of Eq. (16.230 . To calculate the expectation values of Y{f) we try two differnet approaches, 
one where we set 



In — = 

m„ 



(6.24) 



and simply get rid of this term, and a second one where we use the recipe 
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with Inn ~ 0.9 as given in |TY94) . 

The resulting expectation values are shown in Tab. 16.31 for charmonium and in Tab. 16.41 for 
bottomonium. As one can see the size of the corrections is relatively small, but the discrepancies 
between the results with different assumptions are significant. As we do not know a unique way 
to treat the term ln(A/mg) and for the reason that the corrections {Y{f)) are very small with 
assumption Eq. (j6.25p . we consider them to be negligible in our model. 

6.4 Summary 

In this chapter we adopted a fourth-order gluonic quark-antiquark potential by Gupta et al, 
which has been slightly modified by us to guarantee time reversal invariance of the potential. 
This potential is then treated as a perturbative correction. We have predetermined the renor- 
malization scale appearing in the potential based on an argument by Titard and Yndurain. 

In the study of model [CpL*-2] it occurred that the spin-dependent splittings are in general 
too small. For bottomonium the triplet P-state splittings are too small. Remarkably, however, 
this model yields the first time a strong coupling constant < 0.3 for bottomonium. 

In the last section of this chapter we investigated additional terms to the potential by Gupta 
et al, which have been found by Titard and Yndurain. These terms include a fictive gluon mass. 
In the analysis we showed that these additional terms yield only small contributions, but are 
sensitive to assumptions about the fictive gluon mass. Because of this behaviour we decided not 
to take these contributions into account in our models. 



Chapter 7 



Fourth-order Gluonic Potential plus 
Induced Interaction 

7.1 Perturbative Corrections 

In the previous chapters it became clear that on the one hand the induced interaction (Eq. 
(I5.29P ) provides a suitable framework for the triplet P-state splittings. On the other hand 
corrections at 1-loop order to the one-gluon exchange potential (Eq. (|6.1|) ) seem to remedy 
some of the parametrical oddities of the basic model. The next natural step is to merge these 
two models. 

The combination of induced interaction and the fourth-order gluonic potential has an in- 
herent double counting problem. The calculation of Eq. (j6.ip [GR8H ITY94j includes the box 
diagram in Fig. 17.11 One realizes that this boxed diagram is also included in the induced inter- 
action. In principle one would have to subtract this diagram (Fig. 17. ip from the combination of 
Eq. (|5.29p and Eq. (j6.ip . This is not done in the present work, but we correct for this deficiency 
by treating the coupling constant Qj as free parameter. 

7.2 Potential Model 

The potential in this model is composed of the perturbative QCD potential Eq. ()6.ip . the 
induced interaction potential Eq. ()5.29p and a linear confining part: 

^[CpL*pi-2](^-..^) ^ yW(f;p) + V-nd(r;p) + ar . (7.1) 

As for all Coulomb-plus-Linear models the potential used in the Schrodinger equation to calculate 
bound states is 

dry 

Voir) = -^ + CTr, (7.2) 




Figure 7.1: Box diagram included in the compution of Eq. (|6.ip . 
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while the remaining parts of Eq. ()7.ip are treated perturbatively. The mass formula is written 



as: 



+ mP(A;25+1/ .) + mT(A;2^+1;^.) + mP(A;2^+1/,) + M^'^ie'^+H,) , (7.3) 
with the following contributions: 

Cc/,25+1, 33 - 2n/ / In (/iGRr) + 7E 



TT \ \ r 
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(7.4) 
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(7.5) 



The purely momentum dependent terms of Eq. (17. ip are treated according to option 2 of 
Appendix [Dj Employing any other options results in much lower string tensions and higher 
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values for the quark masses, outside the range of empirical values. Figures, tables and quantities 
related to this model are marked as [CpL*pI-2]. 

7.2.1 Charmonium 

The parameters in model [CpL*pI-2] are determined by a fit to the masses of states rjc, J /ip and 
he, and the ratio (Eq. (|3.11|) ). with the additional constraint that the parameters as, 

rric and a are as close as possible to their empirical values. The resulting parameter set is: 

as = 0.295 , 
ai = 1.0 , 

rUc = 1.4056 GeV, (7.6) 
a = 1.04 GeV/fm . 

The model gives values for the string tension a, the strong coupling constant as and the charm 
quark mass nic which are very close to empirical values. 

The reduced radial wavefunctions for the lowest lying bound states are presented in Fig. 17.21 
and Fig. 17. 3^ while in Fig. 17.41 the radial densities for charmonium states are shown. Comparing 
the radial densities of model [CpL*pI-2] with the ones of model [CpL-B], one realizes that they 
are positioned lower in the [CpL*pI-2] plot and therefore the eigenvalues are smaller. This is 
mainly due to the decrease of the string tension by approximately 25 percent. 







Experiment (E"'"04] 


Theory [CpL*pI-2] 


State 
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Mass [MeV] 
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Table 7.1: cc-state masses from experiment and model [CpL*pI-2] (a^ = 0.295, Oj = 1.0, 
a = 1.04 GeV/fm, rric = 1.4056 GeV, M = 3.07 GeV, /iQR = rric, = 3); masses and widths 
are displayed for the experiment. 



54 



Chapter 7. Fourth-order Gluonic Potential plus Induced Interaction 




Figure 7.2: Charmonium S'-state reduced radial wavefunctions calculated within model [CpL*pI- 
2] {as = 0.295, ai = 1.0, a = 1.04 GeV/fm, = 1.4056 GeV, M = 3.07 GeV, hgr = rUc, 
Uf = 3). 




Figure 7.3: Charmonium P- and Z?-state reduced radial wavefunctions calculated within model 
[CpL*pI-2] {as = 0.295, at = 1.0, a = 1.04 GeV/fm, = 1.4056 GeV, M = 3.07 GeV, 
AiGR = '"T'c, = 3). 
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[GeV] [GeV] 




Figure 7.4: Radial densities for charmonium states of model [CpL*pl-2] plotted together with 
used potential (a^ = 0.295, ai = 1.0, a = 1.04 GeV/fm, = 1.4056 GeV, M = 3.07 GeV, 
A*GR = "^o '^Z = 3); base lines of radial densities have been shifted by according E^i. 
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Figure 7.5: Charmonium spectra from experiment and model [CpL*pI-2] (a^ = 0.295, Oj = 1.0, 
£7 = 1.04 GeV/fm, rric = 1.4056 GeV, M = 3.07 GeV, /Uqr = rric, n-f = 3); for experimental data 
masses and widths are shown. 
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Table 7.2: Contributions to the masses of cc-states in model [CpL*pI-2] (og = 0.295, ai = 1.0, 
a = 1.04 GeV/fm, = 1.4056 GeV, M = 3.07 GeV, ^gr = m^, n/ = 3). 



In Tab. 17.11 masses and widths from the experiment are displayed together with masses 
calculated in model [CpL*pI-2]. The experimental spectrum is compared to the one arising 
from model [CpL*pI-2] in Fig. 17. 5[ The data presented in Tab. 17. II shows, that the three masses 
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of our input states r]c, J/tp and he are matched with 1 MeV accuracy. The ratio 

4CpLV-2l(ip)^r ^^(Xc.2)-M(XcaA ^ 5 (7 7) 

VA^(Xc,l)-M(Xc,0)y[CpL*pI-2] 

whicli is used as a further criterion to determine the parameters, agrees quite well with the 
experimental value 

$Exp(ip) ^ ( M(xc,2) - M(xe,in ^ ± 0.005 . (7.8) 

We consider this adequate to constrain the model parameters, as we restricted ourselves to 
certain intervals for three of the parameters and thus the fit is not free. The size of splittings 
of the induced interaction and fourth-order potential did not compensate each other exactly, 
resulting in an enlargement of the triplet P-state splitting by about seven percent. The splitting 
between the singlet and the triplet 25'-states is enlarged, but we match the singlet state rj'^ within 
the uncertainties of the experiment and realize that the triplet state ip' is positioned slightly 
below the L'D-threshold. Therefore the position of ^' could change and the 25'-splitting might 
get a value close to the experimental one. Not surprisingly, the potential model is not sufficient 
to describe the spectrum above threshold. 

In Tab. 17.21 the various contributions to the masses of cc-states are displayed. One notices 
that 

= + Em (7.9) 

leads to values significantly larger than the masses of the respective charmonium states. Hence 
large negative shifts are required. The purely momentum dependent term yields a large 

negative contribution to the masses. An additional large contribution M^'^ for all states comes 
from the radiative correction to the Coulomb part of the potential. The Mg^"^ term of order 
0{m^^) in the interaction also provides large negative mass shifts, but acts mainly on the S- 
states. Together with the contact term M^"^ this influences the string tension parameter in the 
present model as the IS'-states are further separated from the IP-states. 

Tab. 17.21 shows the importance of the induced interaction. The combination Mp + + 
yigifjg 3^22 additional separation between the S'-states and the angular excited states, leading 
to a slightly lower string tension and thus to a value close to the estimate from lattice QCD 
simulations. The induced interaction provides about 40 percent of the singlet-triplet-splitting 
of the S'-states. It also leads to a triplet P-state splitting with the correct structure. This is 
due to the enhancement of the tensor part over the spin-orbit part in the induced interaction as 
compared to the gluon exchange interaction. 
7.2.2 Bottomonium 

The input used to determine the parameters of model [CpL*pI-2] for bottomonium is given by 
the masses of states T(IS') and T(2S'), the center-of-gravity mass C(1P), and the ratio 



$f P(1P) = 2) i Z) [ , (7.10) 



sExp.T _ Mixb2) - M[xbi) 
M{xbi) - M{xm) 

identifying C(1P) with the l^Pi-state of the model. The parameters fullfill the constraint to be 
as close as possible to their empirical values. The resulting parameter set is given by 

as = 0.255 , 

ai = 1.0 , 

rub = 4.768 GeV, (7.11) 

a = 1.08GeV/fm. 
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Figure 7.6: Bottomonium 5-state reduced radial wavefunctions calculated within model 
[CpL*pI-2] {as = 0.255, at = 1.0, a = 1.08 GeV/fm, nib = 4.768 GeV, M = 2mb, hgr = 
Uf = 4). 




Figure 7.7: Bottomonium P- and Z?-state reduced radial wavefunctions calculated within model 
[CpL*pI-2] {as = 0.255, Oi = 1.0, a = 1.08 GeV/fm, rub = 4.768 GeV, M = 2mb, /icR = 2/3™,^, 
Uf = 4). 
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Figure 7.8: Radial densities for bottomonium states of model [CpL*pI-2] plotted together with 
used potential (a^ = 0.255, Oj = 1.0, a = 1.08 GeV/fm, rrib = 4.768 GeV, M = 2mb, hgr = 
2/3mfe, nj = 4); base lines of radial densities have been shifted by according E^i- 
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Figure 7.9: Bottomonium spectra from experiment and model [CpL*pI-2] (as = 0.255, = 1.0, 
a = 1.08 GeV/fm, m{, = 4.768 GeV, M = 2mh, fJ-cR = 2/3mfc, Uf = 4); r]h mass measurement 
by ALEPH collaboration is included separately in the experimental spectrum. 
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Table 7.3: Comparison of experimental masses and widths of 66-states with masses calculated 
with model [CpL*pI-2] (a^ = 0.255, at = 1.0, a = 1.08 GeV/fm, mb = 4.768 GeV, M = 2mb, 
AiGR = 2/3m;,, nj = 4). 



Remarkably, model [CpL*pI-2] implies a value of as for bottomonium which is significantly 
smaller than for charmonium, while the bottom quark mass and the string tension are close to 
their empirical determinations. 

The resulting reduced radial wavefunctions are presented in Fig. 17.61 and Fig. 17.71 while 
in Fig. 17.81 the radial densities of the lowest lying 66-states are displayed. The radial densities 
are drawn at their respective energy eigenvalues to indicate their position in the potential. The 
comparison between Fig. 17.81 and the corresponding figure for model [CpL-B] (Fig. 13. 8p shows 
that the eigenvalues in model [CpL*pI-2] are significantly larger. This is mainly a consequence 
of the smaller coupling and to a lesser degree due to the slight increase of the string tension 
a. 

Experimental masses and widths for bottomonium are displayed together with the masses 
of model [CpL*pI-2] in Tab. 17.31 while Fig. 17.91 presents the spectrum arising from experiments 
together with the one of model [CpL*pI-2]. The value for the mean IP mass is 

C{IP) ^ 9900 MeV . 

The largest deviation for our input states is less than 2 MeV for the T(2S')-state. The ratio 

M{Xbl) - M{Xb0) ) [CpL*pI-2] 
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Table 7.4: Contributions to the masses of 6&-states in model [CpL*pI-2] (a^ = 0.255, aj = 1.0, 
a = 1.08 GeV/fm, mb = 4.768 GeV, M = 2mb, MGR = 2/3mb, n/ = 4). 
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agrees with the experimental value 

^TniP) = = 0.58 ±0.03, 

almost within the uncertainties. Matching <I>(1P) exactly to the experimental value implies that 
the parameters tend to be not as close to their natural values anymore. In particular, the string 
tension increases further after tuning the remaining parameters. Further, the ratio ^b{2P) for 
the 2P-states will be reduced from its value 

calculated with the set Eq. (|7.11|) . which agrees perfectly with the empirical 

= 0.57 ±0.05. 

The contributions to the masses of bottomonium states are displayed in detail in Tab. 17.41 
It is obvious that 

M]:r = '^rnb + Eki, (7.12) 

leads to significantly larger values than the experimental masses. As in the case of charmonium, 
there exist large downward shifts for the bottomonium masses. The purely momentum depen- 
dent term yields large negative contributions. The radiative correction to the Coulomb 
term provides a major negative shift to all states. The remaining terms which provide spin- 
independent splittings, Mg*^^, M^°'^, and M^°'^, act primarily on the S-states and separate 
these states and from the rest. The Mf^^ contribution shows quite a different behaviour for the 
iS-states in bottomonium as compared to charmonium. While in charmonium the magnitude of 
increases with the radial excitation, bottomonium yields the largest contribution for the 
groundstate. This is due to the influence of the Coulomb part of the potential. The induced 
interaction gives rise to about 40 percent of the ^Si — ^Sq splitting. The importance of the in- 
duced interaction for our model manifests itself once more in the structure of the triplet P-state 
splittings. One recognizes that due to the increased weight of the tensor part in the induced 
interaction compared to the gluon exchange interaction, the model agrees with the experimental 
ratios and $b(2P). 
7.2.3 Remarks 

In the effort of describing the quarkonium spectra within potential models, the model [CpL*pI- 
2] represents a substantial improvement compared to the basic model [CpL-B]. It provides a 
detailed description of the spectra below the heavy-light meson thresholds with minimal phe- 
nomenological input and realistic parameter values. Nevertheless some details have to be com- 
mented on. 

There remain questions about the renormalization scale which we predetermined inspired 
by |TY94j . At least for charmonium the link between Os and predetermined scale is not 
obvious. The fact that, for charmonium, as(;U(^) > 0.5, might be a sign that a perturbative 
approach is not justified anymore. 

The induced interaction is only implemented for the IS'-states. In principle one would have 
to implement the induced interaction for all remaining bound states, starting with the induced 
interaction for the IP-states (Tab. 17. 5p . The major argument against considering the exchange 
of other bound states is, that this leads to additional parameters. 

The combination of induced interaction and gluon exchange potential to order O(a^) bears 
the double counting problem mentioned previously. One would in principle have to subtract the 
double counted diagram, which is not yet done in the present calculations. 
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State 
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k^Pi 
k^Po 
k^Pi 
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?/; (7/^75) V' 
? 2kV'(7^(p+p')')V' ? 



Table 7.5: Experimentally determined quantum numbers for the P-states and presumptive 
structures for the t-channel exchange of P-states. 

7.3 Mass-Radius-Relationship 

After studying the spectrum gained via model [CpL*pI-2] we will again test the mass-radius 
relation. In Chapter 13.31 we have already seen that, according to the virial theorem, the energy 
eigenvalues should follow 

„ / 2as 3cr \ ,^ ^ 



resulting in the mass 



Thus one expects that the quarkonium states show an almost linear mass-radius relation. 

To probe the relation of radius and mass one has to decide which masses to take into account 
for the quarkonium states as the perturbatively implemented corrections do not change the 
radius. In Chapter 13.31 we argued that one should take the bare mass 

Mir = '^^<i + Eki, (7.15) 

which does not include any of the perturbatively treated corrections as none of them occur 
explicitely in the virial theorem. In the discussion of model [CpL*pI-2] we realized that the 
values of M^f""^ are about 200 MeV larger than the respective experimental quarkonium state 
masses. A possibility to test the mass-radius relation with masses close to the experimental ones 
is to include the spin-independent corrections of model [CpL*pI-2] into the considered mass 

Mlj = 2mq + Eki + M^" + up + M^""" + M^^ + Mp + Mf°° + M^'^ . (7.16) 

Hence we will investigate the mass-radius relation in two ways, a first with M^^"^"^ as input and 
a second in which M^l is used. 

7.3.1 Charmonium 

It is instructive to examine the radial densities in Fig. 17.41 which demonstrate that the char- 
monium states are dominated by the linear part of the potential and thus a linear mass-radius 
relation emerges. 

In Tab. 17.61 the considered cc-masses and radii within model [CpL*pI-2] are presented. 
Matching straight lines to the data in Tab. 17.61 yields, using M^j^^^^ as input, for the parameters 
of the best fit straight line 

a}"'"' = 2.66 GeV, (7.17) 
1.66 GeV/fm , 



Lbare 
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Figure 7.10: M^f''^ versus (r) for cc-states in model [CpL*pI-2] (a^ = 0.295, a, = 1.0, a = 
1.04 GeV/fm, = 1.4056 GeV, M = 3.07 GeV, jjgr = mc,nf = 3). 
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Table 7.6: Charmonium radii and masses in model [CpL*pI-2] (a^ = 0.295, = 1.0, a = 
1.04 GeV/fm, rric = 1.4056 GeV, M = 3.07 GeV, /icR = mc, n/ = 3). 



while using results in 

Mil = aS^ + 6Si(r), 
a^^ = 2.43 GeV, (7.18) 
b^^ = 1.69GeV/fm. 

The straight line M^^^'^ is displayed together with the input data in Fig. 17.101 The deviations 
from a linear behaviour are small, which seems natural as the viral theorem predicts a linear 
relation between mass and radius in this region. One also recognizes that the value for })°^'^'^ ig^ 
as predicted by the virial theorem, very close to 3/2 cr. 

In Fig. 17.111 the resulting straight line M^ and the input data for the second case are 
presented. One realizes that, while the parameter a^^ is, as expected, about 200 MeV smaller 
than a^^'"'^, the slope parameter b^^ increased slightly comared to the bare mass fit. This is 
understandable if one remembers that some of the perturbatively treated corrections yield an 
additional division between the 5-states and the angular excited states. This additional sepa- 
ration is similar to an increase in string tension and is therfore noticeable via a larger value for 
the slope parameter of the straight line. Further one realizes that the deviations to the linear 
behaviour are considerably larger in this case, but still small enough to justify the assumption 
of a linear mass-radius relation. 
7.3.2 Bottomonium 

For guidance one takes once more a quick look at the radial density plot (Fig. 17. 8p , where one 
can easily behold that the 15-state seems to be special as it is, in contrast to the other states, 
strongly influenced by the Coulomb part of the potential. Consequently a linear dependence of 
mass and radius might only be a good guess for states other than the 15-state. 

In Tab. I7.7l the radii and masses calculated within model [CpL*pI-2] are presented. Matching 
straight lines to the data in Tab. I7.7l vields. using M^f"^^ as input, for the parameters of the best 
fit straight line 

^bare ^ 9.40 GeV , (7.19) 
1.70 GeV/fm , 



Lbare 



while using Mfl results m 



Mil = a^^ + b^\r), 
a^^ = 9.24 GeV, (7.20) 
1.72 GeV/fm , 



LSI 
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Figure 7.12: M^f''*' versus (r) for 66-states in model [CpL*pI-2] (a^ = 0.255, Oi = 1.0, a = 
1.08 GeV/fm, nih = 4.768 GeV, M = 2m;,, /xgr = 2/3m;,, = 4) including the virial theorem 
curve (Eq. ()7.2ip ): IS-state is not included in straight line fit. 




r igure 7.13: M§ versus (r) for 66-states in model [CpL*pI-2] (a^ = 0.255, = 1.0, a = 
1.08 GeV/fm, rub = 4.768 GeV, M = 2mh, /iQR = 2/3mb, nj = 4); IS-state is not included in 
straight line fit. 
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Table 7.7: Bottomonium radii and masses in model [CpL*pI-2] (a^ = 0.255, ai = 1.0, a = 
1.08 GeV/fm, m,, = 4.768 GeV, M = 2mb, ficR = 2/3mfc, n/ = 4). 



where in both fits the IS'-state data has not been included. 

The straight line M^^^^ is displayed together with the input data in Fig. 17.121 In addition 
the prediction of the virial theorem is included in Fig. 17.121 by displaying 

M{{r)) = 2m,-^ + ^. (7.21) 

Fig. 17. 121 illustrates that the deviations of the model state masses from the linear behaviour are 
tiny. Furthermore, even the slope parameter 5^'^''^ is near the prediction of the virial theorem. 
The model data points are lying slightly below the curve inspired by the viral theorem, which 
is due to the fact that one is not able to calculate an expectation value for the hole curve. To 
test this statement we have checked the expectation values of Eq. (j7.2ip for the S'-states and 
found that the values agree with the bare mass values. 

In Fig. 17.131 the resulting straight line M^^ and the input data for the second case are 
presented. One realizes that the deviations from a linear behaviour are significantely larger 
than in the case of bare masses. Compared to the bare mass fit one observes again a slight 
increase of the slope parameter. 

The approximation of linear mass-radius relation for bottomonium is justified for both of 
the discussed cases, as long as the l^-state is excluded. 

7.4 Summary 

In this chapter we have combined the previously discussed models [CpLpI-3] and [CpL*-2] into 
a new model with perturbative corrections based on induced interaction and a fourth-order 
gluonic quark-antiquark potential. The resulting model [CpL*pI-2] incorporates the benefits 
of both models, yielding a quantitatively good spectrum by using reasonable parameters. At 
present this success still comes at the cost of a double counting problem and issues about the 
link between the relevant renormalization scale and the running coupling. 

We found that within model [CpL*pI-2] the strong coupling constants Og for bottomonium 
and charmonium are, for the first time in this work, correlated in the proper way, namely taking 
a significantly smaller value for bottomonium than for charmonium. Furthermore the parameter 
values for the quark masses and the string tension are close to their empirical values. 

The triplet P-state structure of model [CpL*pI-2] coincides, up to some percent, with the 
experimental one for charmonium as well as for bottomonium. 

In a last step we investigated, analogously to the discussion of model [CpL-B], the relation 
between mass and radius. In contrast to the discussion in Chapter we used two distinct 
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assumptions for the considered masses, as the bare masses used in Chapter 13.31 are significantly 
larger in model [CpL*pI-2] than the experimental masses. The analysis for both cases shows 
that a linear mass-radius relation is a good approximation for charmonium in general, and also 
for bottomonium if one excludes the IS-state. 



Chapter 8 

Summary and Conclusion 



8.1 Summary 

In this work, we have developed a potential model which gives satisfactory results for the quarko- 
nium spectra with only four adjustable parameters and minimal phenomenological input. 

Starting out from the conventional Breit interaction combined with a linear confinement 
potential (Chapt. [3]), we argued that the retardation corrections included in the Breit interaction 
are not justified, since the on-mass shell condition is not valid for a bound state. This led us to 
a gg-potential consisting of linear confinement and the one-gluon exchange without retardation 
(Chapt. H]). Such restricted potentials face the following problems: The relative splittings of the 
triplet P-states {Xqjt j = 0,1,2) are grossly overestimated. The string tension for charmonium 
comes out too large. Moreover, the effective coupling Og for bottomonium is significantly larger 
than for charmonium. 

In Chapt. [5] we have presented a possible solution to the problems concerning the triplet 
P-state splitting, by including in the gg-potential contributions from t-channel bound state 
exchanges (induced interaction). Such an additional interaction is a natural consequence of 
crossing symmetry together with the fact that the Schrodinger equation resums ladder-diagrams 
only in the s-channel. However, the problems concerning the string tension and the coupling Os 
remained. 

In the effort to cure also these problems (Chapt. [6|), we made use of a perturbative QCD 
potential by Gupta and Radford, including 1-loop corrections at order 0(0^). We have taken 
special care of time reversal invariance by choosing a more symmetric kinematics in the center 
of mass frame (under time reversal: p — > —p, q ^ q, with q Fourier conjugated to r). This 
procedure led to a reduction of the strength of some contact term by a factor of two. We have 
demonstrated that the inclusion of the fourth-order gluonic potential results in a reduced strong 
coupling as for bottomonium. 

Finally, in Chapt. [7] both types of corrections - induced interaction and two gluon exchange 
- have been combined. Altogether this led to a satisfactory description of the charmonium and 
bottomonium spectra below the heavy-light meson thresholds DD and BB. In particular, such 
an improved potential model is able to reproduce correctly the relative splitting of triplet P- 
states. This comes about through additional interactions of spin-orbit and tensor type. At the 
same time the parameters, quark mass m^, string tension a and strong coupling Og, of the model 
take on values close to those given by QCD determinations and lattice QCD simulations. 

As an application we have investigated the relation between masses and mean radii of the qq- 
states in the model (Chapt. 13.31 fc Chapt. 17. 3p . In each case we find an approximately linear 
mass-radius relation for all states except the 15-state of bottomonium. This underlines that the 
majority of the (/g-bound states is primarily determined by the linear confining potential. 
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8.2 Conclusion 

In this work, we have constructed an improved potential model for heavy quarkonium. Our aim 
has been to get a satisfactory description of the quarkonium spectra with minimal phenomeno- 
logical input. Our corrections beyond linear confinement and one-gluon exchange are physically 
well motivated. They originate either as 1-loop corrections from perturbative QCD or can be 
seen as non-perturbative effects (bound state exchanges in the t-channel) required by crossing 
symmetry. In comparison to other works we did neither allow for an additive constant in the 
confining potential nor did we make any assumptions about its Lorentz-structure (scalar versus 
vector confinement). The conclusions to be drawn from the present work are the following: 

• Potential models with the spin-dependent interaction terms from the one-gluon exchange 
alone are not sufficient to give a quantitative description of the quarkonium spectra (in 
particular fine and hyperfine splittings). Moreover, the parameters - quark mass niq, string 
tension a and strong coupling - in such models take on values incompatible with QCD 
determinations and lattice calculations. 

• Corrections from perturbative QCD such as the fourth-order gluonic potential proportional 
to are essential for an improved description of quarkonium. Their inclusion allows to 
have Os approximately equal for bottomonium and charmonium. 

• The induced interaction (t-channel exchange of gg-bound states) provides a distinct way 
to improve the description of the triplet P-state splittings, different from assuming a 
Lorentz-structure for the phenomenological confinement potential. In the latter case new 
fit parameters are unavoidably introduced. 

• We have demonstrated that a satisfactory description of the quarkonium spectra is pos- 
sible within an improved potential model with reasonable values for the quark mass ruq, 
string tension a and strong coupling Os, and without using additional, uncontrolled phe- 
nomenological input. 

8.3 Outlook 

Even though the potential model constructed in this work provides a satisfactory reproduction 
of the quarkonium spectra, it leaves some open questions. 

The combination of induced interaction and fourth-oder potential potentially includes some 
double counting. However, subtracting the contribution of the particular box diagram (Fig. 17. ip 
cannot be straight forwardly done since it involves a regulating gluon mass. 

The coupling Oi for the induced interaction (15-state t-channel exchange) can be related 
to the asymptotic normalization constant of the wavefunction. This would fix the coupling Oj. 
Furthermore, exchanges of bound states with different quantum numbers (P-states, ...) can also 
contribute to the induced interaction. At present their effect is hidden in the adjusted value of 

A quantitative description of quarkonium near and above the heavy-light meson thresholds 
clearly requires a coupled channel analysis. 

In this work the improved gg-potential has been used exclusively for the description of the 
quarkonium spectra. In the future one should use it also to calculate other properties such as 
decay widths and branching ratios. 



Appendix A 

Conventions and cross sections 



A.l Conventions 

The different conventions employed in quantum field theory only differ in the parameters J^, Ki, 
Li, Ni (i = B, F distinguishes bosons from fermions) [DGH92) . occurring in the normalization 
of fields, 

Jb 

V'(^) = ^l^ms)uip,s)e-'P-- + d\p,s)v{p,s)e'n, (A.l) 
in the normalization of single particle states, 

\k) = W(^)|o), 

\p,s) = LFb^{p,s)\0) , (A.2) 

in momentum space algebraic relations, 

[a{k),a\k')] = KB5^^\k-k'), 
{bip,r),b\p',s)} = KFdrs6^'\p-p'), (A.3) 

and in the normalization of the fermion spinor 

u\p,r)u{p,s) = NF2Ep6rs. (A.4) 

The parameters Jj, Ki, Ni are constrained by the canonical commutation or anticommutation 
relations to obey 

Using the above, the single-particle expectation value of the quantum mechanical probability 
density can be expressed as 

= ^3 i^=B,F). (A.6) 
The parameters of choice used for the calculations in this work are 

Jfj = {27rf/^ J ^ , Kf,j = 1, LFj = {27rf^'y^, Nfj = . (A.7) 
y nij V J J 
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A. 2 Cross sections 

We now want to gain the cross section in a general form, employing the parameters of choice 
defined in App. lA.ll The cross section for two particles to many particles reactions is given by 

da(ai + a2^f) = ^diVfi„ , (A.8) 

where Fine is the flux of incoming particles 

I -> ^1 K\L'^ K2L2 , _ , 

Fine = Pmlvi - V2\ = J^J^l^l - ^2| , (A.9) 

Pfi the transition probability per unit time per unit volume for the process i ^ f 



and dA^fin the density of the final state 

/=! J f 

with Z = Jl^n/j! a statistical factor accounting for the presence of n^t identical particles of type 
k in the final state and X^int ^ sum over internal degrees of freedom. Collecting the pieces, we 
obtain the cross section in the desired form: 

(2-)*^'^'(Ew-E«)E#^.(A-i2) 

/ i int 

where the indices of the sums or products are given by i for particles in the initial state, / for 
particles of the final state and j for particles of initial and final states. 




Appendix B 

Perturbative QCD Potentials 



The short range part of the quark-antiquark potentials in our models is calculated via pertur- 
bative QCD. Therefore, we explicitly calculate one-gluon exchange potentials and dicuss some 
basic problems in their derivation. 



B.l Breit Interaction for Quark- Antiquark Systems 

The Breit interaction has first been deduced by Gregory Breit |Bre291 IBre301 IBre32) for the 
electron-electron scattering process. The according potential includes, apart form a leading 
Coulomb term, relativistic corrections which have their origin in the one-photon exchange process 
and an expansion of the propagator, leading to retardation. Likewise we will rederive the Breit 
interaction for the quark-antiquark scattering process 

qi{PA,SA) +qj{PB,SB) — > QkiPA^s'ji) +qi{PBiSB) , (B.l) 

with four-momenta px, p'x and spins sx, s'^ as indicated, and i,j,k,l = 1,2,3 being color 
indices. In this case the Breit interaction is equivalent to the potential arising from one-gluon 
exchange, including retardation corrections. The scattering process Eq. (|B.1|1 is described, at 
lowest order perturbative QCD, by two Feynman diagrams, a s-channel diagram (Fig. IB.ip and 
a t-channel diagram (Fig. IB.2p . We have to keep in mind, that the purpose of this perturbative 
QCD calculation is to obtain a potential describing quarkonium. As quarkonium states are 
mesons which are color neutral objects, we possess already some knowledge about the color 
wavefunction 

1 ^ 

iMeson) ~ ^ ^ \m) , (B.2) 



V3 , 

m=l 



resulting for Eq. (|B.ip in the conditions 



-^^ij and -^^kl ■ (B.3) 

Taking this condition into account the s-channel diagram does not contribute. Employing Feyn- 
man rules (see e.g. [TWOlj fPSOSj ) the amplitude, corresponding to the s-channel diagram (Fig. 
IB.ip . in Feynman gauge reads 



M 



u{p'a,s'^) {^^97'^) v{p'b,s'b) 



9nu 



(pa + Pb)'^ + ie 



v{pb,sb) ( 'igj"^ ] u{pa,sa) 



^2 



{pA + PB)^ + ie 2 2 



u{p'a, SAh/iviPB, s'b)v{pb, sb)Yu{pa, sa) ■ (B.4) 
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qiiPA,SA 




Figure B.l: One-gluon exchange s-channel diagram for quark-antiquark. 
9fe(PA'«A) \ / 1i{p'b^s'b) 



qi{pA,SA) 




iApb^sb) 



Figure B.2: One-gluon exchange t-channel diagram for quark-antiquark. 



Using Eq. (|B.3p this amplitude vanishes, as 



1 . Af, 



1 



(B.5) 



or in words, the gluon as color-octet does not couple to a color-singlet. This leaves the t-channel 
diagram (Fig. IB.2P to deal with. Employing Feynman rules the t-channel amplitude in Feynman 
gauge is given by 



M 



fi 



u{pa.^a) U91^-f ) u{va.sa) 

9' J 



+ ie 



v{pb,sb) ( 157''^ ) v{p'b^s'b) 



+ ie 2 2 



'^{p'a^ SAhf^u{pA, sa)v{pb,sb)i^v{p'b,s'b) , 



(B.6) 



with = {pa — v'a)^ — {Pb ~ Pb)'^- Remembering Eq. (|B.3p yields the factor 



2 2 -12 

1,1=1 k,l=l a=l 



i,k=l a=l 



a=l 



)'] = 3 • (B.7) 



Thus, for color neutral initial and final states the process Eq. (jB.ip leads to the amplitude 



4 



3 g2 + ie 



u{p'a^ s'AhfMu{pA, sa)v{pb, sb)i^v{p'b, s'^ 



b) 



The computation of the amplitude Eq. (jB.Sp is done partwise. To begin with, we evaluate the 
Dirac currents s'j^'^^u{pA, sa) and v{pb, sb)i^v{jp'q^ s'^), after applying the normalization 
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conventions for Dirac spinors as in Appendix lA.il 



u{pa,sa) 
v{pb,sb) 



' Ea + ruA 
2mA 

'e'^ + tua 
2mA 

'E'^ + mB 
2mB 

I Eb + mB 
2mB 



XSA 

EA+ruA 



XsA 



t ^-p'a 
~^s'j^ E'^+mA 



E'g+mB -^s 

X'Il 



cf 

Xsb EB+rriB ' 



-XsB 



(B.9) 



with energies Ex = -J Px + "f^x ^^'^ quark masses mx ■ Using the auxihary calculations 



Xl/ , X 



s' E' +mA 



CJ^ 
(Ji 



XsA 

_J:PA_ 

EA+ruA'^^A 



t S-pk t 
Xs'^E'^+niA^''' Xs>/t 



c\ S-pB ^4 \ / Ui 

XsBEB+m.B'' XsB J \ (J 



XsA 

_i:PA_ 

EA+ruA ■^^A 
/ S-pk c 



E'g+mB ^s'g 



t / S-p'a 
'^^'a [e' +mA 



■ai + cr, 



■ PA 

Ea + ruA 



Xs 



•A 1 



(B.IO) 



X%cri, X. 



cf a-pB 
SB EB+niB 



O'-PB ,,C 

E'g+niB^s'g 

xl> 



/ysB 



<y -PB 

Eb + rnB 



CTj + (Ti 



E' 



'Xl 



B 



mB 



{a ■ e){a ■ x) + {a ■];){&■ e) 



€■ X + ia ■ {ex x) + y ■ e + if? • (y x e*) 
e- [y + x) + \e- a X (y - x) , 



(B.ll) 



one can easily verify 



u{p'A,s'Ah^u{pA,SA) 



v{pb,sb)i^v{pb,Sb) 



u\p'a^s'j^)u{pa,sa) 



u\p'A^s'j^)-iQ^u{pA,SA 



2mA 



'E'^ + mA IEa + niA 



2mA 




' E\ + mA IEa + ruA 



2mA 



X' 



2mA 

x' 



1 I i^-pKii^-PA) \ ^ 

^ ^ {E'^+m)A{EA+mA) J XsA 



E'j^+rnA ^ " EA+mA ) XsA 



7t <^-PA 



( 1 _i_ {o-pA{o-pa) i ^ 

'a\ {E'^+mA){EA+mA) J XsA 



PA 



Ea+itia 



+ 



PA 



hniA 



+ i(7 X 



PA 



PA 



E'.+mA EA+mA 



SA 



2m B 



E'j^ + mB IEb + tub 



2m B 



(B.12) 
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for the Dhac currents. For the further calculation we substitute 



PA = Pi + 2^ 
1 

2^ 



PB = P2 



PA=pi - ^q, 
PB=P2 + i^q, 



SA = Si , 
SB = S2 , 



mA = mi, 

VflB = 1712 ■ 



(B.13) 



We choose these kinematics to get a coordinate space potential which is invariant under the 
usual time reversal transformation 



t 

f 

Pi 
J 



-Pi , 
-J. 



(B.14) 



The kinematics-choice is essential due to the fact that, while the amplitude is time reversal 
invariant as long as it is expressed in Lorentz-covariants, its nonrelativistic expansion and the 
gained potential are not (see Appendix IB.3j) . With new variables the timelike and spacelike 
parts of Eq. (fRT2]) read 



u\p'a:s'a)u{pa,sa) 



vHpb,sb)v{p'b,s'b) 



EA + mi E'^ + mi | 



2mi 



2mi 



Ea + nil 
2mi 

Eb + m2 
2m2 



2mi ^'^'i 



Ea + mi Ej+mi | 



2mi 



X' 



2mi ^^'^ 
^ E'B + ^\ .ci 

2m2 



1 + 



1 + 



{a-{pi-\q)){^-{pi + \q)) 
{EA+mi){E'j^ + mi) 

{pi - \q) ■ jpi + y) 

{EA + mi)iE'^+mi) 
ia- {{pi - y) X {pi + 



{EA + mi){E'^ + mi) 



1 + 



Pi 



{EA + mi){E'j^ + mi) 



Xs 



Xsi 



Xsi 



\ P2 - - '^^ ■ {q ^ P2) \ c 

{EB + m2){E'^ + m2) i 



'^{p'a^s'a)iu{pa,sa) 



v{pb,sb)iv{p'b,s'b) 



Ea + mi E'^ + mi 
Y 2mi 

Pi + y 



X 



1 



2mi 

t 



+ 



Pi 



2mi '^^'i 

Xs2 [■ 



2m' 



EA + mi ' E'^ + mi 
2pi -iaxq] Xs^ + O 

32 - i(T X q] xl> + O 




Pi + jq 
Ea + mi 



(B.15) 



For the spacelike part of the currents we have neglected corrections of order 0{m^ ). These 

irrlpr (Km^^m,^ 



and 0{m^ "^m - ^) contributions in 



order 0{m- ^) terms in the currents give order 0{m^ "^m- ^ 

the amplitude. These terms are negligible as we expand the amplitude to order 0{m'^'^m^'^). 
The spin of particles is given by the expectation values 



Sl 5s^s{ 



xl'^Xs^ 



(B.16) 
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of the Pauli matrix a. In the case of antiparticles, we have to deal with charge conjugated PauH 
spinors 

S2K,s', = X%<yx\- (B.17) 

Thus, we have to relate the spin to charge conjugated spinors. The charge conjugation of 
wavefunctions is given by 

C = 17^7°. (B.18) 
Employing this we can express x'^^^X'^ with Pauli spinors 

{x'^^x") = xV2a<72X* = -(x^^x)* (B.19) 

where 



(72<T(72 = (72 \ = -B* , (B.20) 




is utilized in the last step. Knowing that the spin is a real quantity, we recognize that 

(X'^^X') = - (x^<?x) , (B.21) 

and therefore 

S2 <5,,,. = - (xL^X^) • (B.22) 
Setting the pieces together, we obtain for the amplitude 

^ {EA + mi){E'^ + m^) j ^'^^^^ ^ ^ {Eb + m2){E'j, + m^) ) 
= -N-A A Pi -ig^-2igi-(gxpi) p| - Iq"^ + 2\S2 ■ (g xp2) \ 

t / 2pW^2il^ , ( 2P2 -id^q\ e , J: 



2. (_4p'-^ . + 4isi • {q X P2) - 4is2 • {q x pi) 

4mim2 

1 



-4r • .?2) + 4 (^-- si) (g • S2) + O 



Ea + itii E'.+mi IEb + m2 EL+m2 



2mi y 2mi V 2m2 y 2m2 
4mf Xmf 
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where we have expanded the energies in the respective last step. In the fohowing the 5g ,,' are 
suppressed. To find the potential corresponding to Ai/i, we have to relate the non-relativistic 
cross section from Schrodinger theory with the cross section computed in Appendix IA.2[ We 
recall that for the non-relativistic case, the differential cross section is related with the Born 
scattering amplitude via 



(27r)3(27r)2 



/ 



\fBl 

\Vrel\ 



(B.24) 



Applying our conventions to the cross section given in Appendix lA. 21 we end up with 



dan = dVid^p^ 



2 2 

m|m2 



Pi 



f 



\Vrel\ 



(B.25) 



Knowing that the Born scattering amplitude is equal to the Fourier transform of the potential, 

d^r 



/b 



-ip -r 



(2^)3 

we gain the relation between potential and amplitude 



17111712 



Expanding the energies the relation reads 



M 



fi 



1 



-Viq;pi,P2) 



± 1 



Pi 



^2 



P2 



Imf Smf 2m\ 8m| 



2 / 1 

+ o' 



M 



(B.26) 



(B.27) 



(B.28) 



The sign has to be determined by the fact that the interaction is repulsive in the case of the 
electron-electron one-photon exchange, giving a minus sign. Combining Eq. ()B.23P and Eq. 
(Ib:28]) leads up to order 0{ ^171- ^) to 



(2vr)^ 



:V{q;pi,P2) 



45^ 



1 + 



_ 2[si ■ {q X pi) 



4mf 
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_^ 2is2 • (q x P2) 



4m| 
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4mim2 



3(g2 + ie) ^ 

-4pip2 + 4isi • {q X P2) - 4is2 • {q x pi) - Aq^ (si • S2) 4 (q- si) 



(B.29) 



S2, 



Since we want to obtain a nonrelativistic potential we have to expand the 4-vector in the prop- 
agator, 
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q'' 




= (E'a- 


- EaKEb - E'j,) = 
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1 {pi ■ q){P2 ■ q) 






q'^ 77117712 



(#1 • q){p2 ■ q) 



m\m2 
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+ 0( — 



1 

777'^ 



(B.30) 



resulting in 
1 



(2^)= 



1 + 



-\q^ - 2isi • (g x px) 



^ 4 (pi • q){p2 ■ q) ^ 

(p 4m 17772 



+ 
1 



\ra\7a2 



-\q^ + 2is2 • (g X P2) 

4777^ 



-4pi^2 + 4isi • (g X ^2) 



-4is2 • {q X pi) - 4g2 (si • 5*2) + 4 (g • si) (g • 5*2; 



(B.31) 
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The Fourier transformation back into coordinate space 

V{f;pi,P2) = / d^qe~"^'''V{q;pip2) , 



(B.32) 



of the various terms in Eq. (|B.3ip is given by 
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W3 
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-iq-r 



[ r 
1 



d^qqiqje-''^''^ = -ViV,U 
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(B.33) 
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47rr2 ' 



W" 



AW = -5(2) (^) , 



-— , [/" = , AC/ = --^ 



(B.34) 



and the computation of ViVjW being performed explicitly in Appendix O Collecting all the 
pieces we obtain the Breit interaction for quark-antiquark scattering 



+ - 



_4a,^ 2^,^(3) 
3r 3 
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2a., 



3mim2 



Pi •P2 (f •pi)(f -^2) 
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3mim2 
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^5(3)(,-.)(,-'^ . + 3(.l-f)(.-2-f)-.-'i-.-'2 



(B.35) 



3^3 



o o + ((r X ^li) • S2 - (r X |)2) • Si) 



rn-im2 



B.2 One-Gluon Exchange without Retardation 

The calculation of the Breit interaction which is equivalent to the one-gluon exchange with retar- 
dation corrections, and the one-gluon exchange without retardation differ only in the treatment 
of the propagator. Therefore, we may as well start the calculation of the one-gluon exchange 
without retardation on basis of Eq. (|B.29|) . Instead of Eq. ()B.30p we use 

1 ^ 1 9o 




'4 ' 

(B.36) 
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for the propagator, leading to the momentum space potential 



(2vr)^ 
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(B.37) 



+ -r-^ ( - 4piP2 + 4isi • {q x P2) - 4is2 • {q x pi) - 4^^ (si • 5*2) + 4 {q- si) {q- 5*2) 

4mim2 V 

Employing the Fourier transformation of the various terms of Eq. given by Eq. ()B.33p . 

yields for the potential in coordinate space 



V^'\r;pi,P2) = + '-^6^'\n ( 4 + 4 ) + , 

or 6 V w-T "^2/ omim2 r 



1 1 



2 ' 2 
m\ m 



Aus pi ■ P2 



4a., 



2>m\m2 
3r3 
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(B.38) 



o o + {{r X Pi) ■S2-{rx P2) • Si) 



mim2 



The label y*^^^ is chosen as the diagrams taken into account to compute the potential possess 
two vetices. 



B.3 Choice of Kinematics 

The kinematics-choice is essential in the derivation of potentials from perturbative QCD due to 
the fact that, while the amplitude is invariant under the usual time reversal transformation as 
long as it is expressed in Lorentz-covariants, its nonrelativistic expansion and the gained potential 
are not. Starting from Eq. ()B.12p we want to show the importance of the kinematics-choice 
for the potential and the effect of retardation, by explicit recalculation of a quark-antiquark 
interaction with different kinematics. 



Non time reversal invariant kinematics-choice 

For this case we choose the kinematics: 

PA=Pi, PA=Pi-q, SA = si, mA = mi, 

Pb=P2, PB =P2+q, SB = S2, TUB =m2 . (B.39) 

For the timelike and spacelike parts this leads to the expressions 
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(B.41) 
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(B.42) 



Employing these expressions in the amplitude Eq. ()B.8p . remembering Eq. (1B.16|1 and Eq. 
(inil]), gives 
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The relation between potential and amplitude 
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and we obtain 
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(B.46) 



Using Eq. (|B.36p for the expansion of the propagator, the potential in momentum space reads 



(2vr) 



1 V 
Viq;pi,P2) = -i7^ 



+ 



1 + 
1 



-^q"^ -2isi ■ {qx pi) -ig^ + 2is2 • x ^2) 



4m| 



+ 



4m2 



4mim2 



-4pip2 - 2pi ■ q + 2p2 ■ q + q + 4isi • x p2 



4is2 • (g X pi) - 49"= (si • S2) + 4 (g • si) (g • 5*2; 



. (B.47) 



As one can see there are spin-independent terms of structure pi ■ q. These structures will after 
Fourier transformation result in terms which are not invariant under the usual time rever- 
sal transformation. The parameterization yields an additional term assigned with a factor 
{mim2)~^ . 



Time reversal invariant parameterization 

For the case, which leads to a time reversal invariant potential, we substitute the variables in 
the amplitude by 



PA = Pi + , PA = Pi - , SA = si , mA = mi , 
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PB = P2 + ^q , SB = S2, mB = m2. 



(B.48) 



As this is exactly the choice used in in App. IB.ll and lB.21 we skip the details and perform only 
the last steps. Employing Eq. ()B.36P for the propagator 
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-4pip2 + 4isi • (g X P2) - 4is2 • (g x pi) - 4g (si • S2) + 4 (g • si) (g • S2 



One recognizes that there are no terms of the structure described in the previous paragraph in 
Eq. (IB30]) . 
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The effect of retardation 

The previous paragraphs discussed the influence of the kinematics-choice on potentials without 
retardation corrections. Taking into account retardation in the same way as for the Breit 
interaction (App. IB.ip . the influence of the kinematics-choice vanishes as both parameterizations 
result in 
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This is a consequence of the distinct expansions for the energies in 
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For the non time reversal invariant case we obtain, 
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while for the time reversal invariant parameterization it eventually reads. 
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Though the statement, that retardation compensates the influence of the kinematics-choice, is 
true for the One-Gluon exchange, it does not hold for the Lorentz-vector exchange of massive 
particles. 



Appendix C 

Computing V^VjVI^(r 



As one has to proceed with caution in computing an expression for 'Vi'VjW^r), depending 
reclusively on derivations of W with respect to r, we wih show the calculation in detail. First, 
we search for the most general ansatz for two derivatives with two indices, and eventually find 

ViVjW = (a6ij + b^) -W + (c6ij + d^) W" + (e5ij + /^) AVF . (C.l) 



The powers of r in the denominator are chosen in this way to make sure that the coefficients 
a, . . . , / are dimensionless. Distinguishing the case of i = j from i 7^ j we multiply our ansatz 
by using that 

b.jhj = tr(%) = 3 , (C.2) 

to arrive at 

AW = (3a + b) + (3c + d) W" + (3e + /) AW . (C.3) 

This yields the following conditions for the coefficients 

b d 1 f 

a = --, c=--, e = ---, (C.4) 

which, introduced into our ansatz, result in 

ViVjW = (^b^W + dW" + fAW^ {^-^ - ^5i^ + ^5ijAW . (C.5) 

The remaining constraint needed to identify our coefficients are gained by investigation of the 
case i 7^ j. Remembering that 

ViW(r) = -W'(r) , (C.6) 
r 

we find 

V.V.^ = 5^ (w" - Iw') , i^j. (C.7) 



i^j. (C.8) 



\ r 

Then again, using the ansatz Eq. (ICSh yields 



ViVjW = ^ {b\w + dW" + /AVF^ 



Through comparison of the last two equations we eliminate the remaining unknowns, giving: 
V,Vj^(r) 



W"ir) - -W'ir) 
r 



^ - \^^^ + \h^W{r) . (C.9) 
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Expectation Values of Momentum 
Dependent Terms 

The computation of quarkonium masses within the model includes perturbatively treated terms. 
One class of terms are purely momentum dependent terms which are given by a Hamiltonian of 
structure 

^MD ^ /(r)(^ap2 + 6(f-p)(f -p)^ , (D.l) 

with an arbitrary function /(r) and a, b being constants. For the evaluation we have to translate 
the Hamilton function into the corresponding Hamilton operator. As the operators for p and 
r do not commute the transition from Hamilton function to Hamilton operator is not unique, 
meaning that there is no unique ordering for the operators corresponding to p and r. Hence we 
have to look at several orderings of operators, with the prerequisite that the resulting expression 
has to be hermitean. 



D.l Options for {H^^) 

In this work we choose one of the three options presented in the following for the momentum 
dependent terms. 



Option 1 

In the first option the ordering is as such that the operator corresponding to the function /(r) 
in Eq. ()D.1|1 is either completely left or right of the p operators. Therefore, the expectation 
value of Eq. (|D.ip in option 1 is given by a superposition of 

{fir)f)i = -\jd\{i,*V^f{T)^ + rf{rW^) 
= -\j ((VV)/WV' + V'*/(r)VV) 
= - j d\V*/(r)VV, (D-2) 
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and 



{f{r){f ■p){r -p))^ 



d-^r ij* f{r){f ■ V)(f • V) + (f • V)(f • V)/(r) ip 



dV 
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In terms of the reduced radial wavefunction u{r) these expression transform to 

dV f{r)i)*V'^i) 



(D.3) 



(/(r)p')i 



(u(r))^ 



+ , 



{f{r){r ■p){r ■p))i 



d?' /(j') u{r)u"{r) 
/ dV /(r)V'*^^ 



(D.4) 



(D.5) 



Option 2 

For the second option we write down a more symmetric form, where the operator correspond- 
ing to the function /(r) in Eq. (jP.ip stands in between the p operators. The expectation value 
of Eq. (|D.ip is then a superposition of the expressions 



{f{r)p')2 



dV V*V/(r)VV' = / d^r /(r) I V T , 



{f{r){f ■ p){f ■ p))2 = - j <rr ^* f -V f{r)f -V i) = j (Tr f{r)\f -V ^\ 
The equivalent expressions in terms of the reduced radial wavefunctions are given by 



(D.6) 
(D.7) 



{f{r)p 



dV/(r) I V V 
dV/(r) I 

dV /(r) I 

dr r'^f{r) 

dr f{r 



d^r f{r) 



V 





2 




+ 




2 




+ 



n(r) 
r 

u(r) 



u{r) 



r xVYi. 



u'{r) n(r)^ ^ ^ n(r) 



dl7 



(D.8) 



{f{r){r ■ p){r ■ p))2 



dV/(r) 
d^r /(r) 



r • V V' 

5 / u(r) 



dV/(r) 
2 



9_ 

dr 



dr \ r 



dr /(r) [u (r)) - 2 ^ — 



(D.9) 



D.2. Momentum Dependent Terms of the Breit Interaction 
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Option 3 

In addition to tiiese two basic options we also employ an option motivated by the Weyl- 
prescription [BCP94| 



1 



{/(OpV}h/ = 7{{/(0,p'},p^} • 



W 4 

This option is essentially the arithmetic mean of options 1 & 2. 



(D.IO) 



Although in principle any other normed mixing of options 1 &: 2 would be possible, we pick, 
for the momentum dependent terms within our model, one of the three options above. 



D.2 Momentum Dependent Terms of the Breit Interaction 



The momentum dependent term of the Breit interaction (Eq. IB.35P is a unique case as option 
1 and option 2 and therefore any superposition of them are equivalent. In the following this 
feature of the Breit interaction is shown. 

The momentum dependent term of the Breit interaction is given by 



Tj-MD 
-"Breit 



2as 
3m^ 



(D.ll) 



The expectation value of Eq. (jP.lip in options 1 and 2 in terms of the reduced radial wavefunc- 
tion (see App. ID.ip read 



A 



2as 



{p^ + {f-p){r-p)) ^ 



^ / drU2u{r)u"{r) 



2u{r)u'{r) {u{r)) 



+ 



2-/(/ + i; 



B 



2as 



(p^ + (f ■ p){r -p)) 



2as 
2>m? 



J dr^ (^-2{u'{r)f + 4 



u'{r)u{r) {u{r)y 



2 + /(/ + l)H . (D.12) 



Comparing the two options one finds that they differ only by a total derivation: 

v2 



A-B 



2a., ^dri (4l^-6 "(^^;^("^ 2(n-(r))V2.(rK(r) 



2>m? 

4as 



dr 



d_ d u{r) 

dr \ dr r 



4as 
2>m? 



(r))^ u{r)u'{r) 



.(D.13) 



r=0 



The boundary condition which the radial wavefunctions must fulfill at r = and the utilization 
of a power series ansatz for u{r) near the origin give (see e.g. |Fri06j ) 



n(0) 



0, 



iimnir) = lim 

r-»0 r^O r 



(D.14) 



Furthermore, we know that in our models the linear part of the potential is dominant for large 
r. Thus, the reduced radial wavefunctions features for r — > oo the asymptotic behaviour of 
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solutions for a linear potential which are given by Airy-functions, yielding the conditions 

lim u{r) = , 

r^oo 

lim u'{r) = . (D.15) 

r— >oo 

Considering these properties of the reduced radial wavefunction we can verify 

'v?{r) u{r)u'{ry 



0, 



(D.16) 



r=0 



implying that options 1 and 2 are equivalent. 
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